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CHAPTER  I  -  INTRODUCTION 

I.  A  GENERAL  BACKGROUND 

A  composite  solid  propellant  consists  of  finely  divided  oxidant 
and  metal  particles  uniformly  dispersed  in  a  matrix  of  polymer  binder. 
The  most  commonly  used  oxidant  is  ammonium  perchlorate  (AP)  and 
some  of  the  fuel  binders  used  are  polyurethane*  polybutadiene-acrylic 
acid,  polysulphide,  polyester,  etc.  In  the  present  investigation  only 
non  metalized  propellants  with  AP  particles  of  a  single  size  (unimodal) 
were  considered. 

In  the  past  few  decades,  a  considerable  number  of  experimental 
and  theoretical  investigations  have  been  directed  towards  the  study  of 
the  steady-state  burning  of  composite  propellants.  The  purpose  of 
these  investigations  has  been  to  or  edict  or  to  allow  willful  modifica¬ 
tion  of  the  performance  of  AP  composite  propellants  in  order  to  avoid 
the  time  consuming,  repetitive,  and  costly  experimental  tests.  In 
spite  of  these  efforts,  the  present  models  are  imperfect  and  are  not 
capable  of  correctly  representing  the  burning  mechanism  of  composite 
solid  propellants. 


2 


1.  B  BURNING  MECHANISM  CF  A  COMPOSITE  PROPELLANT 

The  burning  mechanism  of  a  composite  solid  propellant  is  very 
complex  due  to  the  heterogeneous  composition  of  the  propellant. 
Several  spatial  dimensions  several  phases,  and  several  simulta¬ 
neous  and/or  successive  physical  and  chemical  processes  are 
involved. 

During  burning,  gases  evolve  from  the  propellant  surface  and 
mix  and  react  exothermally  in  the  gas  phase.  The  heat  generated  by 
the  reactions  is  transferred  to  the  surface  and  produces  further 
vaporization  of  the  surface  material.  Previous  studies  (1,  2,  3,  4) 
have  indicated  that  the  following  physical  and  chemical  processes 
are  important  in  the  burning  of  a  composite  solid  propellant  : 

(i)  heat  conduction, 

(ii)  surface  pyrolysis  of  the  fuel  (endothermic), 

(iii)  surface  dissociative  sublimation  of  AP  into  ammonia  and 
perchloric  acid  (endothermic), 

(iv)  vapor  diffusion, 

(v)  AP  decomposition  flame  (exothermic  reaction  between 
ammonia  and  perchloric  acid),  and 

(v)  final  flame  (exothermic  reaction  between  pyrolysed  fuel  and 
AP  decomposition  products). 

..  .  -  


Figure  1-1  illustrates  the  burning  mechanism  schematically. 

The  linear  rate  at  which  the  surface  regresses  macroscopierlly 
is  a  characteristic  of  the  burning  and  is  designated  as  the  burning 
rate  of  the  propellant  (r).  The  desired  end  results  of  steady-state 
mechanistic  models  are  the  prediction  of  r  and  of  the  influences  on  r 
by  pressure  (P),  oxidant  particle  size  (doy),  initial  temperature  (TQ) 
and  weight  fraction  of  oxidant  (mixture  ratio)  (w). 

I. C  PURPOSE  AND  SCOPE  OF  THE  PRESENT  INVESTIGATION 

I.  C.  1  P^JP^se_ 

The  purpose  of  the  present  investigation  is  to  develop  an 
improved  steady-state-burning  model  for  AP  composite  solid 
propellants.  The  need  for  such  a  model  based  on  reasonable 
mechanistic  grounds  arises  from  the  inadequacy  of  the  existing 
models  to  describe  the  experimentally  observed  effects  of  all  the 
important  variables  on  the  burning  rate. 

I.  C.  2  _Scop_e 

The  presently  reported  work  consists  of  : 

a.  review  of  prior  experimental  and  theoretical  investigations 
of  the  steady-state  burning  of  AP-eomposite-solid  propellants, 

b.  new  experimental  and  theoretical  investigations  concerned 


FINAL  FLAME 


SCHEMATIC  OF  THE  BURNING  PROCESS  OF  A  OTPOSITE  SOLID  PROPELLANT 


with  mixture-ratio  influences  on  burning  rates, 

c.  assessment  of  existing  comprehensive  models  by  means  of 
numerical  diagnostic  tests,  and  finally, 

d.  formulation  of  a  modified  granular  diffusion  flame  model, 
based  on  reasonable  mechanistic  grounds  and  testing  of  the  model 
against  experimental  data. 
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CHAPTER  II  -  REVIEW  OF  PRIOR  INVESTIGATIONS 

In  the  present  chapter,  prior  experimental  and  theoretical 
investigations  of  the  steady- state  burning  of  AP  composite  solid 
propellants  and  their  analogs  are  reviewed. 

n. A  EXPERIMENTAL  INVESTIGATIONS 

A  great  deal  of  effort  has  been  devoted  in  the  past  to  study 

burning  rate  behaviour  of  AP  composite  propellants  with  respect  to 

pressure  (P),  oxidant  particle  size  ( d  ),  initial  temperature  (T  ), 

ox  0 

and  mixture-ratio  (w).  Both  actual  propellants  and  their  analogs 
have  been  used  in  these  investigations. 


II. A.  1  investigations  with_  Propellants 


A  detailed  review  on  this  subject  can  be  found  in  references  3 


o 


5  H  C'ViC  f4  f'T'Jv  °  hr' o  i » monxr  nf  v  ooi  1 1  f  c  Tro  (Titron  V»ol  mu 

w  «  WX  XV/i  U  V4  XXllltilX  J  VX  w»»V  .»  vx;  v><.  vw  v  i  w»*  wa  w  <•  • 


II.  A.  l.a  Pressure 

In  general,  r  is  found  to  increase  as  P  increases.  Based  on  the 
variations  of  r  with  P,  three  pressure  domains  have  been  defined 
as  follows  : 

(i)  low  pressure  domain  (below  10  atm. ), 

(ii)  moderate  pressure  domain  (10  to  100  atm. ),  and 

(iii)  high  pressure  domain  (above  100  atm. ). 
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These  demarcations  are  not  rigid  as  they  deform  as  other  parameters 
are  varied.  These  different  domains  are  attributed  to  the  domination 
of  the  overall  burning  process  by  different  phenomena  (e.  g,  chemical 
reaction  and  diffusional  mixing  at  different  pressures). 

II.  A.  l.b  Oxidant  Particle  Size 

Bastress  (6)  has  done  a  systematic  study  of  the  influences  of 

dQX  on  r  and  found  that,  in  general,  r  increased  with  decreasing  dQX. 

The  influence  of  d  on  r  has  also  been  found  to  be  coupled  with 
ox  r 

variations  in  other  parameters,  such  as,  oxidant  weight  fraction  (3). 

II.  A.  1.  c  Initial  Temperature 

The  burning  rate  usually  increases  as  the  initial  temperature  is 
increased.  A  considerable  interest  is  shown  at  present  in  under¬ 
standing  the  influence  of  initial  temperature  on  the  burning  rate  (7). 

II. A.  l:d  Mivhire  Ratio 

A  measure  of  the  mixture  ratio  in  this  study  is  taken  to  be  the 
weight  fraction  of  oxidant  (w).  The  influence  of  w  on  r  has  been 
virtually  ignored  by  the  previous  investigators.  Only  very  limited 
data  is  available  with  regard  to  the  mixture  ratio  influences  and  no 
attempt  has  been  made  previously  to  interpret  these  data.  Therefore, 
in  the  present  investigation,  the  available  data  have  been  plotted  as 


,  w... . . . 


r  vs.  w  (none  of  the  previous  studies  shows  a  plot  of  this  type). 
Figures  II-l,  2,  and  3  show  plots  of  r  vs.  w  for  various  pressures, 
particle  sizes,  and  different  types  of  fuel. 

The  following  observations  can  be  made  from  these  figures  : 

(i)  The  burning  rate  (r)  increases  monotonicslly  with  the 

mixture  ratio  (w)  (all  the  w's  considered  are  less  than  the  stoichio 

# 

metrically  correct  w,  w  ^). 

(ii)  The  mixture-ratio  sensitivitv  of  the  burning  rate,  , 

aw 

varies  with  P,  dQX,  w,  and  the  type  of  fuel  used. 

(iii)  r  vs.  v;  curves  in  figuresll-l  and  2  suggest  a  tendency  to 
attain  a  maximum  r  as  w  increases,  though  the  available  data  are 
not  extensive  to  ascertain  this  trend. 

II.  A .  2  Jnvesti gations  with_Analo gs 

Analogs  of  AP  propellants  were  used  in  the  past  to  study 
burning-rate  variations  with  changing  P,  dQX,  T0,  and  v;  (4,  9,10). 
The  advantages  cited  for  using  analogs  rather  than  true  propellants 
are  that  the  analogs  are  easy  to  fabricate  and  that  they  allow  the 
composition  to  be  varied  at  vail,  which  is  not  possible  with  actual 
propellants. 

? - 

ws(.  is  defined  here  as  that  w  which  assumes  the  products  of 

combustion  to  be  HgO,  CO2,  SO2,  HC1,  and  N2. 


12 


The  three  types  of  analogs  used  in  the  past  are  : 

(i)  Pressed-samples  (4,9), 

(ii)  Loose- granule  burners  (10),  and 

(iii)  Porus-plug  burners  (10). 

In  the  present  study,  additional  loose-granule  tests  were  used  to 
investigate  mixture-ratio  infii  ences. 

II.  A.  2.a  Parametric  Influences 

The  results  obtained  from  analog  studies  show  the  same  influenc- 
f  )f  pressure,  oxidant  particle  size,  and  initial  temperature  on  r 
as  those  shown  by  actual  propellants  (4,  9, 10).  However,  the  analog 
studies  provide  additional  information  with  regard  to  the  mixture- 
ratio  influences  on  r.  Unlike  propellants,  analogs  can  be  fabricated 
and  tested  over  a  wide  range  of  mixture  ratios.  Figures  If-4,  5, 
and  6  show  r  vs.  w  curves  from  previous  studies  (4,  10).  These 
curves  are  useful  in  studying  the  influences  of  fuel,  oxidant  particle 
size  and  pressure  on  the  trends  of  r  vs.  w  curves  (  i,  e.  the  shifts  in 
the  mixture  ratio  that  correspond  to  maximum  burning  rate, 

Wr(max)* 

Figure  II-4  shows  that  for  mixtures  containing  AP  and  poly¬ 
styrene  (PS),  the  mixture-ratio  yielding  maximum  burning  rate, 
wr(max)>  *s  011  ^e  ^el  ric^  s^e  °£  the  stoichiometrically  correct 
mixture  ratio  (w  t),  whereas  for  mixtures  containing  AP  and  pciy- 


BURNING  RATE  (CM/SEC) 


FIGURE  I M  :  IfFUUNCE  OF  FUEL  ON  THE  POSITION  OF  .TOMJM 
BURNING  RATE  (LOOSE-GRANULE  BURNERS) 


MIXTURE 


5H 

140-320  f* 


P  *  10  ATM 


methylmetliacrylate  (FMM),  wr^max'  1S  near  wst*  Hence,  it  is 
concluded  that  the  type  of  fuel  used  in  the  propellant  has  an  influence 
on  the  position  of  the  maximum  burning  rate. 

The  influence  of  the  oxidant  particle  size  on  the  position  of 
maximum  burning  rate  is  shown  in  Figure  II- 5.  For  a  larger  oxidant 
particle  size,  ^‘r(rnax)  is  found  to  be  on  the  fuel  rich  side,  whereas 
for  a  smaller  particle  size,  wr (max)  is  found  to  be  near  the  stoichio- 
metricallv  correct  mixture  ratio,  w 

st 

The  influence  of  pressure  on  the  position  of  maximum  burning 

rate  is  shown  in  Figure  II-6.  At  low  pressures  wr(max)  is  less  thuu 

w  and  as  the  pressure  is  increased  w  .  shifts  towards  w  .. 
st  r(max)  st 

The  experimental  data,  shown  in  Figure  IT-6,  were  not  taken 
beyond  w  and  hence  the  existence  of  the  maximum  burning  rates 
'with  w  variations  at  high  pressures  could  not  be  conclusively 
confirmed.  Part  of  the  present  investigation  involves  measuring 
regression  rate  for  a  wide  range  of  variations  in  w  at  different 
pressures  using  loose- granule  burners.  The  details  and  results  of 
such  measurements  are  given  in  Chapter  111. 


II.  A.  S  Classification^ of  Propellants 

AP  composite  propellants  have  been  classified  by  Steinz  (2)  into 
two  categories:  those  that  do  and  those  that  do  not  allow  correlation 


)7 


of  r  vs.  P  data  by  means  of  Summerfield's  equation  (1),  viz. 

I  _  a  .  b  ...  II -1 

T  '  v  T'y 

where  a  -  chemical -time  parameter,  and 

b  -  diffusion-time  parameter. 

Both  a  and  b  are  assumed  to  be  independent  of  P. 

The  agreement  (or  lack  of  agreement)  between  the  r  vs.  P  data 

P  2/3 

and  Equation  II-l  can  be  determined  by  plotting  the  data  as  ^  vs.  P 
and  checking  whether  such  plots  result  in  a  straight  line.  This 
method  of  classification  is  illustrated  in  Figures  II-7  and  8. 

Figure  II-7  shows  examples  of  propellants  which  do  correlate  via 
Equation  II-l,  and  Figure  II-8  shows  an  example  of  a  propellant 
which  does  not  follow  Equation  II-l.  The  corresponding  r  vs.  P 
curves  are  shown  in  Figures  II-9  and  10  respectively.  It  can  be  seen 
clearly  from  these  figures  that  Summerfield's  equation  is  a  good 
representation  for  those  propellants  with  burning  rates  (r)  increas¬ 
ing  monotonically  with  P  and  with  the  slope  of  r  vs.  P  curves 
decreasing  steadily. 

In  the  present  study,  only  those  propellants  for  which  data  can 
be  correlated  with  Equation  II-l,  are  considered. 


AP/LP5  PROPELLANT 


PROPELLANTS  WHICH  OBEY  EQUATION  II-l 


BURNING  RAIB  (CM/SEC) 
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n.  A.  4  ^J]x^rimenta.l_  Investigations 

Based  on  a  review  of  prior  experimental  investigations*  the 
following  observations  can  be  made: 

(i)  In  general,  there  is  a  lack  of  comprehensive  experimental 
data  available  which  deal  with  mixture-ratio  influences*  In  order  to 
assess  existing  and/or  new  models  of  propellant  burning  more 
comprehensive  data  is  needed. 

(ii).  None  of  the  propellants  which  follow  equation  n-1  were 
tested  previously  with  more  than  two  variations  in  the  mixture  ratio 
(2).  Since  a  number  of  propellants  obey  Summerfield's  equation, 
additional  r  vs.  w  data  should  be  obtained  for  such  propellants. 

H.  B  THEORETICAL  INVESTIGATIONS 

Theoretical  models  describing  the  steady-state  burning  of 
composite  solid  propellants  can  be  classified  into  two  categories 
as  follows: 

(i)  simple,  analytical  diagnostic  aids,  and 

(ii)  comprehensive  models. 

The  simple,  analytical  diagnostic  aids  neglect  many  of  the 
physical  processes  and  chemical  reactions  for  the  purpose  of 
mathematically  tractable  solutions.  They  lay  their  emphasis  on 
specific  processes.  The  following  models  belong  to  this  category: 


(i)  Phalanx  flame  model  (11), 

(ii)  Premixed,  laminar  flame  model  (12)^ 

(iii)  Sub-surface  reaction  model  (13), 
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(iv)  Thermal  layer  theory  (14), 

(v)  Columnar  diffusion  flame  model  (15),  and 

(vi)  Quasi— laminar,  diffusion,  and  monopropellant  flame 

theory  (16), 

Since  these  models  neglect  many  of  the  presumably  important 
processes  and  do  not  explain  the  burning-rate  beheviour  completely, 
they  are  not  considered  further  in  the  present  investigation. 

The  comprehensive  models,  on  the  other  hand,  are  based  on  more 
complete,  detailed  mechanistic  grounds.  They  allow  prediction  of 
the  burning-rate  behaviour  with  respect  to  variations  in  important 
parameters. 

Hence,  in  the  present  study,  the  following  comprehensive  models 
have  been  considered; 

(i)  Granular  diffussion  flame  models  (GDF), 

a.  Original  GDF  model  (1,  6). 

b.  GDF  model  with  collapsed  AP  decomposition  flame 
(GDFC)  (2). 

c.  GDF  model  with  distended  AP  decomposition  flame 


(GDFD)  (2). 


— ’T-*;sgy  ^  swwfSW 
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(ii)  Heterogeneous  reaction  model  (HR)  (17, 18), 

(iii)  Multiple  flame  model  (MF)  (19,  20). 

II .  B .  1  Granular  D iffus ion_Flam e_ Model_  (GD F) 

II.  B.  l.a  Description 

In  the  general  formulation  of  the  GDF  model  (2),  the  gasification 
of  the  solid  surface  was  considered  to  be  driven  by  the  conductive 
heat  feed  back  from  a  two-stage  flame  in  the  gas  phase. 

The  first-stage  flame  was  considered  to  be  the  AP  decompostion 
flame,  and  the  second-stage  flame  was  considered  to  be  the  result 
of  the  reaction  between  pockets  of  fuel  (granules)  and  the  oxidant 
rich  products  of  AP  flame.  It  should  be  noted  that  there  is  no 
experimental  evidence  for  the  existence  of  the  fuel  pockets,  and  this 
assumption  has  been  seriously  questioned  (3). 

A  schematic  description  of  the  GDF  mode]  is  shown  in 
Figure  li-ii. 

n.  B.  1  .b  Assumptions 

Itemized,  the  major  assumptions  of  the  GDF  model  are  : 

(i)  One-dimensionality  (planar  surface  and  reaction  zones  in 
the  gas  phase). 

(ii)  Steady-state. 
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FIGURE:  II-ll  :  DESCRIPTION  OF  TOE  GDF  MODEL 
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(iii)  Specific  physico-chemical  processes: 

The  physical  processes  considered  are  heat  conduction  and 
granular  diffusion. 

The  chemical  reactions  considered  are  fuel  pyrolysis, 
dissociative  sublimation  of  AP  into  ammonia  and  perchloric 
acid,  AP  decomposition  flame  (first-stage  flame  resulting 
from  a  second  order,  premixed  reaction  between  ammonia 
and  perchloric  acid),  and  granular  diffusion  flame  (second- 
stage  flame  resulting  from  the  diffusional  mixing  and  reac¬ 
tion  between  fuel  pockets  and  oxidant  rich  products  of  AP 
flame). 

(iv)  Linear  temperature  distribution  in  the  gas  phase. 

(v)  Uniform  thermal  properties. 

II.B.Ijc  Equations 

T?r»r\vY\  fV»  o  Hftcr*rinfinri  fhr*  oocnmnfinno  r.ohHh  in  fVin  frwtvnilo  — 

A.  «  V  AAA  wiiv  V  A  W  t  J  v  A  A  V  A  V_/  A  A  A  A  Vt  V*  *  W  kV  W  A  4  A  ^  »*V  lAku<  AA  4AA  V4  V  AAA  Vi.  V  A  V/ A  .  4  A  »  A  44 

tion  of  the  GDF  model,  the  following  equations  can  be  written  (2). 

An  energy  balance  at  the  surface,  for  the  first -stage  flame  gives: 

V^-)*  PPrlc3CTs-TjtwaAP  +  ci-w)a,] 

*  •  .  .  *U‘" 


....... ..  ... . 
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An  energy  balance  at  the  first-stage  flame,  for  the  second-stage 
flame  gives  : 

ypr  [c, C5 -t„ )  +  w qap o-w) - w a hap] 

II- 3 


where 


Ag  -  thermal  conductivity  of  gases, 

Cg  -  specific  heats  of  solid  and  gases, 

Tq  -  initial  temperature, 

Tg  -  surface  temperature, 

-  first-stage  flame  temperature, 

T2  -  second-stage  flame  temperature, 

-  first-stage  flame  height, 

L0  -  second- stage  flame  height, 

J  -  propellant  density, 

r  -  burning  rate, 

w  -  weight  traction  of  oxidant, 

Q^P  -  heat  of  sublimation  of  AP, 

Q  -  heat  of  pyrolysis  of  fuel,  and 

r 

AH  -  heat  of  reaction  of  AP  flame. 

AP 


An  Arrhenius  expression  is  assumed  to  describe  the  surface 
pyrolysis  of  the  propellant  and  is  given  by: 
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r 

-  As  (“Es/rtj) 

...  II-4 

where 

A 

s 

pre-exponential  factor, 

E 

s 

activation  energy,  and 

R 

universal  gas  constant. 

II .  B .  2  Or i  gir* al  _GDF  _Mod  el 

II .  B .  2 .  a  Additional  Assam p t ions 

In  the  original  CDF  model  (1),  the  following  additional 
assumptions  were  made  : 

(i)  Constant  surface  temperature,  T  . 

s 

(ii)  An  AP  decomposition  flame  (first-stage  flame)  occuring 
entirely  ("collapsed”)  at  the  regressing  surface. 

II.B.2.b  Formulation  of  thp  Model 

The  collapsed  AP  flame  assumption  reduces  Equation  II-2  and 
Equation  11-3  to  a  single  equation  : 

X9  (5ul)=  ?? r  [c9 Os -t.)  +waAT+o-w)«f-wAHAJ 

...  n-5 

In  Equation  II-5,  L£  is  unknown  and  is  expressed  by  considering  the 
detailed  structure  of  the  granular  diffusion  flame.  To  approximate 
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Lg,  two  extremes  of  low  and  high  pressures  were  considered.  At 
low  pressures,  chemical  reaction  was  considered  to  be  controlling; 
at  high,  pressures,  diffusional  mixing  was  considered  to  be  controll¬ 
ing. 


Letting  L2  ^  denote  the  granular  diffusion  flame  height,  when 

chemical  reactior  is  controlling  (at  low  pressures),  the  chemical 

mass  conversion  rate  in  a  zone  of  thickness  L„  .  may  be  written 

2,  ch 

as: 


?„  re„  =  L B  ?„  exf>  (- Ex/k t)  ...  n-6 


where 


ch 


B 


pre- exponential  factor,  and 


E  -  activation  energy  for  the  gas-phase  reaction 

Li 

Substituting  the  value  for  Lg  Cj1  from  Equation  II-6  in  Equation 
II-5,  an  expression  for  r  ^  has  been  obtained  (2)  : 

I  c  cr  _  f  +  w  +  n-  W)  Q ,  -  WAH  r, 

-  y  I  3  -  -  •  't\r  '  TT 

trt"  P  4  A9  C^-Tj) 


*(^)  exKEy2RTJ 

where  oc  is  a  constant. 

On  the  other  hand,  considering  high  pressures  and  control  of 
burning  by  the  diffusional  mixing  of  the  fuel  pockets,  an  expression 
for  Lg  has  been  obtained  (2)  : 


-  . . - 


L; 


ft 


By  substituting  the  value  for  in  Equation  II- 5,  the  burning 


rate  at  high  pressures,  r  ,  has  been  obtained  (2)  as: 


i  |  c3  cr$  -t0)  +•  w  Qap  +  o-^,)  Of  -  ^  ^  ^ap 

few 3  » 


x  TkVt  \ 

V  -pv3  TjV®  / 


Vr*~7o> 


II-9 


With  the  two  asymptotic  forms  for  the  burning  rate,  given  by 
Equations  II-7  and  9,  the  burning  rate  at  an  intermediate  pressure 
was  expressed  by  a  relation  which  reduces  *0  one  of  these  limiting 
forms  at  either  extreme  of  pressure. 


The  following  form  was  proposed  in  reference  (1): 


I 


I  l 

7.  *  c 


rr  1  a 

...  XX  X  v. 


c* 


dl*f 


or,  substituting  from  Equations  IT-7  and  11-9  : 


=  % 


C9  (Ts  ~T»)  +  w  GW  +■ 


^9 


r 


X  | exp  ( ^/KrJ  *  ^5-— *5 


11-11 


■:-.-*S**ti*IU*H**- . 


31 


or 


and 


...  n-12 


•py» 


where 

0-  =  ff 


c9  ^rT|>)+ wQap  4(i-w)QrwAHw 


*^9  C  *i"  %  ) 
X  £c<T2  e*P  (E*/rT*)J 


...  11-13 


b  -  £, 


Cg  CTS  -Ta )  4-W  QAp  4*  (l-W}  Qj  ~  WaH 
AgC-Ti-TJ) 

*4 


Af» 


f  PC4ox)t2*  ? 

1  V*  1 


11-14 


II.  B.  2.  c  Discussion 


In  the  preceding  equations,  the  burning  rate  appears  in  a 
transcendental  form  requiring  an  iterative  solution.  The  principal 
equations  involved  are  Equations  II-4  and  II - 1 1 .  The  two  unknowns 
are  r  and  T,  the  burning  rate  and  the  surface  temperature,  respec- 

b 

tively. 

In  the  original  formulation  of  the  CDF  model,  these  equations 
were  not  solved  numerically.  The  assumption  that  the  surface 
temperature,  T  is  constant  allowed  the  Equation  II-4  to  be 

k) 

neglected,  and  the  burning  rate  was  related  to  the  pressure  by 
means  of  Equation  II  - 12  only;  this  equation  is  sometimes  known  as 
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Summerficld  equation.  The  validity  of  the  GDF  model  was  said  to  be 
established  by  showing  tha  t  experimental  burning-rate  data  (r  vs.  P) 
could  be  correlated  by  the  Summer  field  equation  for  many  propellants 
(1,2,6). 

II.  B.  3  GDFC  _Model 

Steinz  (2)  has  modified  the  GDF  model  by  allowing  the  surface 
temperature  to  vary  with  pressure,  in  accordance  with  Equation  II-4, 
but  retained  the  collapsed  AP  decomposition  flame  assumption.  This 
model  was  termed  the  "granular  diffusion  flame  model  with 
collapsed  AP  decomposition  flame"  (GDFC)  model. 

II.  B.  3,  a  Resulting  Equations 

The  equations  of  the  GDFC  model  are  those  of  the  GDF  model: 

r  =  A s  C~  )  ...  xi -4 

JL-  Q  -  CS  ~T0^  +  w'c^y  t  0-W)  q7-  w  A  Hat. 

r"  A3(7i-T5) 


e*b/£r.'|  +  P(-U 


Equations  II-4  and  11-11  contain  two  unknowns,  r  and  T  .  They 

s 

can  be  solved  numerically  by  using  the  meihod  described  below. 


II.  B.  3.  b  Method  of  Solution 


The  input  values  required  to  solve  Equations  II-4  and  II- 1 1  are 
classified  as  r 


(i)  environmental  input  values  (P  and  TQ), 

(ii)  compositional  input  values  (w,  dox,  ,  T2,  QAp,  Qf,  AHAp), 

(iii)  kinetic  input  values  (E  ,  E{ ,  etc. , ),  and 

(iv)  empirical  input  values  (A  ,  c(,  and  ). 

The  empirical  input  values  are  obtained  as  follows  : 

,  P  2/3 

(i)  Determine  values  for  a  and  b  from  a  plot  of  —  vs.  P 

(ii)  Choose  a  pressure  (matching  point),  determine  the  corres¬ 
ponding  value  for  r,  and  select  an  appropriate  value  for  T  . 

(iii)  Use  Equations  II-4  to  determine  A  . 

s 

(iv)  UseEquations  II- 1 3  and  11-14  to  determine  ©(and  (J. 

Equations  II-4  and  II- 1 1  can  now  be  solved  numerically  for  various 
other  pressures. 


TT*i  »  v»  a  TT  1  O  aaaviIKaa  fU  a  w  aI-'U  a/4  it  a  Arl  I* a!  trn  fVt  f~\  T?0  *r\  ol 

X'  VC-jV-AA  'C.  n“JLW  UUOUi  IWV-O  WIW  iuuk.auu  WO'-U  tv  OV1V  v  W*v  t-U-/ J.  mvv*v-i 


equations. 


II.  B.  3.  c  Results  and  Discussion 

In  Reference  2,  numerical  solutions  were  obtained  for  a  typical 
propellant  by  considering  only  pressure  variations.  The  predicted 
r  vs.  P  results  were  shown  to  have  the  same  form  as  that  of 
Summerfield  Equation  (Equation  11-12)  in  1.  0  to  100  atm.  range. 


FOR  GDF  MODEL 


,y  -•  j*  Jtliira''l-T  i^^i  Tri  ^> 


The  predicted  results  were  not  compared  with  any  experimental  data. 
The  matching  point  needed  to  obtain  the  empirical  constants  A„,  o(  , 

b 

and  ^  ,  was  chosen  arbitrarily  at  a  low  pressure.  The  possible 
effects  of  alternative  selections  for  the  matching  point  (at  different 
pressures)  was  ignored.  Furthermore,  the  influences  of  oxidant 
particle  size,  mixture  ratio  and  initial  temperature  on  the  burning 
rate  were  not  discussed. 

Therefore,  as  part  of  the  present  investigation,  the  following 
comprehensive  numerical  tests  were  made  to  assess  the  GDFC 
model  : 

(i)  A  comparison  of  the  predicted  GDFC  results  with  experimen¬ 
tal  results  (rvs.P). 

(ii)  An  investigation  of  the  effect  of  the  matching  point  selected. 

(iii)  An  investigation  of  the  influence  of  mixture  ratio  on  the 
burning  rate  as  predicted  by  the  GDFC  model. 


II.  B.  4  GDFD  Model 


The  collapsed -AP-flame  assumption  made  in  the  GDFC  model 
was  found  by  Steinz  (2)  to  be  invalid  at  low  pressures.  Hence, 

Steinz  proposed  a  GDF  model  with  "distended"  AP flame  (GDFD).  The 


GDFD  model  is  supposed  to  be  valid  at  low  combustion  pressures. 
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II.  B.  4.  a  Resulting  Equations 

As  a  consequence  of  the  above-mentioned  assumption,  the  GDFD 
model  is  described  by  Equations  II-2,  II-3,  and  II-4.  In  order  to 
solve  these  equations,  the  flame  height  of  the  first-stage  flame 
must  be  specified. 

Considering  mass  conservation,  can  be  written  as  (2)  : 

,  r  ...  n-15 

l  “  p~~  *  la  p 

J9 


where  is  the  reaction  time  of  the  AP  flame.  Assuming  the 

AP  flame  to  be  a  premixed  and  second-order  reaction,  7"'  was 

AP 

found  to  be  (2)  : 


6-5  x  »o~6 

? 


(sec) 


...  n-ie 


The  flame  height,  L  ,  has  been  taken  to  be  the  same  as  that 

Cl 

- 1  .1 — l  it.  • 

tuiwiuw  cu  in  tuts  kju r  v  muutjji. 


The  following  equations  resulted  after  simplification  (2)  : 

f  =  As  exj*  (-E5/RTf) 


l_  S.p+('-W)Qfj 


n-i7 
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,  _  f  X3  CVt;  ) 

Lifp  {C3^-r»)+w  «t4P  +<i«w)Qf -W4Ha, 

...  11-18 

II.  B.  4.  b  Method  of  Solution 

The  resulting  equations,  viz.  Equations  II-4,  17,  and  18  of  the 

GDFD  model  contain  three  unknowns  r,  T  ,  and  T..  The  numerical 

s  1 

procedure  used  to  solve  these  three  equations  is  the  same  as  that 
described  in  Section  II.  B.  3.  b  and  also  in  Figure  n-12. 

II.  B„  4.  c  Results  and  Discussion 

The  results  reported  in  Reference  2  for  the  GDFD  model  were 
obtained  for  the  same  typical  propellant  used  for  the  GDFC  model 
calculations  described  earlier  (Section  II.  B.  3.  c).  The  burning 
rates  predicted  by  the  GDFD  model  in  the  1.  0  to  100.  0  atm.  range 
were  found  to  be  less  than  those  predicted  by  the  GDFC  model. 

Thp  GDFD  model  rpsults  wprr-  not  Compared  with  any  pYpprinipptnl 
data.  The  influences  of  matching  point  selection,  oxicia.it  particle 
size,  mixture  ratio,  and  initial  temperature  were  not  reported 
in  Reference  2. 

Therefore,  as  part  of  the  present  investigation,  the  following 


numerical  tests  were  made  to  assess  the  GDFD  model: 
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(i)  A  comparison  of  the  predicted  GDFD  results  with  experi¬ 
mental  data  (r  vs,  P). 

(ii)  An  investigation  of  the  influences  of  matching  points. 

(iii)  An  investigation  of  the  influence  of  mixture  ratio  on  the 
burning  rate  as  predicted  by  the  GDFD  model. 

These  tests  and  their  results  are  discussed  in  Chapter  IV. 

II.  B .  5  H eter ogeneouj _Rea.ction_  Model_  JHR) 

II.  B.5.a  Description  of  the  HR  Model 

Hermance  (17, 18)  developed  a  model  of  composite  -  solid  - 
propellant  burning  by  incorporating  effects  of  propellant  surface 
heterogeneity  and  of  heat  released  at  the  surface  due  to  a  heteroge¬ 
neous  reaction  between  the  oxidant  pyrolysis  products  and  solid 
fuel  binder.  This  model  can  be  used  to  predict  known  combustion 
characteristics  like  the  burning  rate,  surface  temperature,  etc,, 
a  priori  ,  using  only  the  propellant  composition  and  component 
decomposition  characterist.es  as  input  data. 

Figure  TI-13  describes  schematically  the  structure  of  the  surface 
and  the  sites  of  heat  generation  considered  in  the  HR  model. 
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ILB.  5.  b  Major  Assumptions 


(i) 

Geometry  : 

solid  phase 

one  dimensional 

surface 

three  dimensional 

gas  phase 

one  dimensional 

(ii) 

Steady  state. 

(iii)  Uniform  thermal  properties. 

(iv)  Specific  physico-chemical  processes  : 

The  physical  processes  considered  are  heat  conduction,  and 
diffusion . 

The  chemical  reactions  considered  are  fuel  pyrolysis, 
dissociative  sublimation  of  AP  into  ammonia  and  perchloric 
acid,  exothermic  reaction  between  fuel  binder  and  products 
of  AP  sublimation  (heterogeneous  reaction),  and  exothermic 
reaction  between  final  AP  decomposition  products  and  fuel 
vapors. 

II.  B.  5.  c  Resulting  Equations 

Due  to  the  complexity  of  the  model  and  the  numerous  equations 
involved,  a  detailed  review  is  not  given  here.  Hence,  only  the  final 


equations  are  listed  below  : 


mw*\ 


r  =  «,  ^(-i^)  + 

X  e*\>  (  “  ft>*  ”  ^st 


b,p 


es 


L 


zi_ 

\r 


t1*! 


E  «, 


) 


n-19 


e„  = 


b*  +  (^)  [^-^1  ^['^sr1'] 

L  2 

+  fhr  ex{3  (-  g|5)t  a3  e*|>(-(3f  &*  e  f(£)  ] 


...  11-20 


JL 

04  '-0S  +a3  [*-  Cxl»EPflf  **({)*]} 


11-21 


where 


and  the  constants  ap 


and  other  parameters  involved  are  defined  in  References  17  and  18. 


W*5W5,!li!|«'!P 
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II.  B.  5.  d  Method  of  Solution 

The  above  equations  contain  three  unknowns,  viz.  r,  Tq,  and  Tp 
Numerical  solutions  were  obtained  by  an  iterative  procedure 
described  in  Reference  17. 

II.  B„  5.  e  Results  and  Discussion 

The  results  of  the  HR  model  were  compared  with  experimental 

data  of  AP/LP3  (6)  propellants  only.  The  model  has  been  tested 

comprehensively,  and  qualitative  agreement  has  been  found  between 

the  predictions  of  the  model  and  experimental  data  with  regard  to 

the  influences  of  P,  d  ,  w  and  T  on  r  (17,  lb).  Since  the  experi- 

ox  o 

mental  data  available  for  mixture  ratio  variation  were  confined  to  a 
narrow  range  ol'  w  (0.  6  to  0.  7),  the  results  of  tic  HR  model  were 
not  considered  beyond  wsO.  7.  Hence,  an  additional  test  useful  for 
evaluating  this  model  would  be  to  obtain  more  results  beyond  w:  0.7 
and  to  study  the  predicted  trends  of  r  vs.  w  curves. 

One  of  the  most  severely  criticised  assumptions  of  the  HR  model 
is  the  heterogeneous  reaction  considered  to  be  occur ing  in  the  cusps. 
There  is  no  experimental  evidence  for  the  occurrence  of  such  a 


reaction. 
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II.  B.  6  Multiple  Flame  Model  JMF) 

II.  B.  6,  a  Description  of  the  MF  Model 

The  multiple  flame  model  has  been  based  on  a  flame  structure 
surrounding  individual  oxidant  particles.  Three  sepei  ate  flames 
have  been  considered  to  supply  the  heat  required  for  surface 
gasification.  They  are  as  follows  : 

(i)  A  primary  flame  (PF)  surrounding  the  AP  particle  caused 
by  the  reaction  between  a  fraction  of  AP  sublimation  products 
and  a  fraction  of  fuel  vapors.  Both  kinetics  and  diffusion 
are  assumed  to  be  controlling. 

(ii)  The  AP  decomposition  flame  (AP  flame)  controlled  by 
kinetics  alone  . 

(iii)  A  final  flame  (FF)  caused  by  the  reaction  between  AP 
decomposition  products  and  fuel  vapors  and  is  controlled 
by  diffusion  alone. 

In  Figure  11-14,  a  schematic  description  of  the  MF  model  is 
given. 
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II.B.S.b 


Major  Assumptions 


(i)  Geometry  : 

Solid  phase  -  one  dimensional 

Surface  -  three  dimensional 

Gas  phase  -  primary  and  final  flames-  three 

dimensional,  and  AP  flame  -  one 
dimensional. 

(ii)  Steady  state. 

(iii)  Uniform  thermal  properties. 

(iv)  Specific  physico-chemical  processes  : 

The  physical  processes  considered  are  heat  conduction  and 
diffusion. 

The  chemical  reactions  considered  are  fuel  pyrolysis,  AP 
sublimation,  primary  flame,  AP  decomposition  flame,  and 
final  flame. 


H.B,6.c 


Equat  ions 


Mass  conservation  at  the  surface  gives 


K*  ?Pr=*C(|f)+  *♦(&)  ••• n-22 
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where 


mass  flux  of  the  propellant, 
mass  flux  of  the  oxidant, 
mass  flux  of  the  fuel, 


ox  surface  area  ratio  of  oxidant  to 

So' 

total  propellant,  and 
S 

f  surface  area  ratio  of  fuel  to  total 

so 

propellant 

n^ox  and  m^  have  been  expressed  in  terms  of  pre- exponential  factors, 

A  and  A  ,  and  activation  energies,  E  and  E  .  and  the  surface 
ox  f ’  b  ’  ox  f 

temperature,  T  ,  as  given  below  : 
s 


*0,  =  Aox  exf>  (-E°*/4?T5)  •••  n-23 


rtif  =  At  e*t>  (~  E+/-RTs  ) 


...  11-24 


s  s 

ox  and  _£  were  expressed  as  functions  of  w  and  other  factors 
S0  So 

which  influenced  the  surface  geometry. 
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An  energy  balance  at  the  surface  was  given  by  : 


Ts  =  T0  _  w  |t  -  (t-w)  Si  +  c^jwf  %?  e-|>  (-  f  J 


L  ^ 


+  ^£  &*M'S*F)]  +  fe  ^£f  e*(>(-f*f) 


...  n-25 


where  Q’s  represent  the  heat  of  vaporization  of  solids  and  heats 
released  by  various  flames,  pFis  the  fraction  of  oxidan*  reactants 


*i 


that  react  in  the  primary  flame  (PF),  and  l[  $  are  the  nondimensional 

*  .*» 

flame  stand-off  distances  of  various  flames,  (T=  rEiHxj. 

*  3  A 

The  flame  heights,  x  ’s,  were  obtained  for  kinetic -controlled 

flames  and  diffusion-controlled  flames  in  the  following  forms  : 

.  « 


m 


rs 


...  11-26 


(kinetic-controlled  flame)  ^  ^ 
where  k  is  the  pseudo-rate  constant  and  Sis  the  reaction  order 


x  =  f(^',  o . )  11-27 

(diffusion-controlled  flame) 

1 

where  is  the  nondimensional  flame  height  given  by  Burke  and 


Schumann's  thin-flame  analysis  (22), 
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ILB.60d  Method  of  Solution 

The  equations  of  the  MF  model  contain  basically  two  unknowns, 

r  and  T  ,and  can  be  solved  numerically  by  an  iterative  procedure, 
s 

The  burning  rate  behaviour  with  pressure,  oxidant  particle  size, 
mixture  ratio,  and  initial  temperature  has  been  studied  (19,  20). 

II.  B„6,  e  Results  and  Discussion 

The  formulation  of  the  MF  model  is  found  in  two  publications 
(19,  20).  The  major  assumptions,  the  resulting  equations,  and  the 
method  of  solutions  were  the  same  in  both  these  references,  but  the 
input  values  used  were  different. 

The  following  comments  are  made  with  regard  to  these  two 
presentations  : 

(i)  Reference  19  contains  a  complete  set  of  input  data  needed 
for  numerical  calculations  of  the  MF  model;  Reference  20 
contains  incomplete  input  data. 

(ii)  The  results  reported  in  Reference  19  are  poorly  supported 
by  experimental  data.  The  results  of  Reference  20  are 
somewhat  improved. 

In  order  to  assess  this  model,  a  complete  set  of  input  data  used 
in  Reference  20  has  been  obtained  from  one  of  the  authors(23)}  and 
further  tests  have  been  made  on  the  MF  model  as  part  of  the  present 


4S 


study. 

Mixture  ratio  effects  considered  in  References  19  and  20  were 
confined  to  a  narrow  range  of  w  (0. 6to0. 7  ),  In  order  to  study 
the  trends  of  r  vs.  w  predicted  by  the  MF  model  more  calculations 
are  needed,  and  new  numerical  tests  made  on  the  MF  model  in  the 
present  investigation  are  discussed  in  Chapter  IV. 

II,  C  CONCLUDING  REMARKS 

The  preceeding  review  of  the  prior  experimental  and  theoretical 
investigations  of  the  burning  of  AP  propellants  has  established  the 
need  for  futher  consideration  of  the  influences  of  mixture  ratio,  and 
the  need  for  further  assessment  of  the  existing  models  by  more 
numerical  tests. 

In  Chapter  III,  new  experimental  and  theoretical  investigations 
undertaken  in  the  present  investigation  with  regard  to  the  mixture- 
ratio  influences  are  presented. 

In  Chapter  IV,  an  assessment  of  the  important  theoretical  models 
is  made  based  on  a  number  of  new  numerical  test. 
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CHAPTER  III  -  PRESENT  INVESTIGATIONS  INVOLVING 

MIXTURE-RATIO  INFLUENCES 


HI.  A  BURNING  RATE  MEASUREMENTS  ON  PROPELLANT 
ANALOGS 


The  influences  of  fuel  and  oxidant  particle  size  on  the  mixture 
ratio  for  maximum  burning  rates,  wr^ma  were  studied  previously 
by  Bakhman  et  aL  (4),  and  McAievy  et  al.  (10)  by  varying  w  over  a 
wide  range  (beyond  w  on  the  oxidant  rich  side).  But  the  data  obtain- 

SI 

ed  by  Bakhman  et  al.  (4)  to  study  similar  influences  of  pressure, 
were  not  taken  beyond  w  .  The  occurrence  of  maxima  in  r  vs.  w 
curves  has  also  not  been  well  confirmed  above  10  atm,  (see  Figure  • 
II-6).  Hence,  burning  rates  have  been  measured  in  the  present 
investigation  over  a  wide  range  of  w  (beyond  wgj.),  for  different 
pressures,  by  using  a  loose- granule  burner.  The  purpose  of  these 
measurements  is  to  supplement  the  previously  obtained  data  (4)  and 
to  verify  those  results. 

The  fabrication  of  the  loose-granule  burner,  the  apparatus  used, 


and  the  techniques  used  to  measure  the  burning  rate  are  described 
in  detail  in  Reference  10,  and  hence,  only  new  experimental  results 
are  discussed  at  present. 
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III.  A.  1  Results  and  Discussion 


The  results  of  the  present  investigation  together  with  those  of 

reference  9  are  plotted  in  I  gure  III- 1.  It  is  seen  clearly  that 

v  .  shifts  from  the  fuel-rich  side  towards  tne  stoichiometrically 
r(max) 

correct  mixture-ratio  as  pressure  is  increased.  This  behaviour  is 
in  conformity  with  the  results  of  Bakhman  et  al.  (4). 


III.  B  EXTENSION  OF  A  SIMPLIFIED  PREMIXED  FLAME 
MODEL  TO  ENCOMP ASS  COMPOSITE  PROPELLANT 
BURNING 


In  the  past,  solid  propellant  burning  models  were  based  on 
simple  extensions  of  laminar,  premixed  flame  theories  to  predict 
burning  rate  characteristics,  such  as,  the  dependence  of  the  burning 
rate  on  pressure,  flame  temperature,  etc. ,  (3, 12,  16).  Since  the 
mixture  ratio  influences  the  structure  of  a  premixed  flame  jia 
flame  temperature  a  simplified,  premixed,  and  laminar  flame 
model  has  been  formulated  as  part  of  the  present  investigation  to 
explore  mixture-ratio  influences  on  the  burning  rate  of  composite 
solid  propellant. 


BURNING  RATS  (CM/3EC) 


FIGURE  III-l  • :  INFLUENCE  OF  PRESSURE  OH  THE  POSITION  OF  PWXIUM 
BURNING  RATE 
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III.  B.  1  Description  of  tbe  _Model 

Figure  m  -2  illustrates  the  simple,  premixed,  laminar -flame 
model  considered.  A  single  premixed  flame  in  the  gas  phase  is 
considered  to  release  heat  which,  conducted  to  the  solid  surface, 
causes  regression  of  the  solid  surface.  The  temperature  of  this 
flame  is  assumed  to  be  equal  to  the  adiabatic  flame  temperature 
of  the  propellant. 

III.  B.  2  Major_  Assumptions 

(i)  The  solid  surface  and  the  flame  arc  planar  and  parallel  to 
each  other  (one  dimensional). 

(ii)  Steady  state. 

(iii)  Uniform  thermal  properties. 

(iv)  An  Arrhenius  type  expression  describes  the  surface 
pyrolysis  of  both  AP  and  fuel. 

(v)  A  premixed,  laminar  flame  of  second  order  exists  in  the 
gas  phase, 

(vi)  The  temperature  distribution  in  the  gas  phase  is  linear, 

ill.  B.  3  Equations^ 

The  surface  pyrolysis  of  the  propellant  is  assumed  to  be 
described  by  an  Arrhenius  type  expression  given  below: 


FIGURE  1 1 1-2  :  THE  SIMPLIFIED,  PREFIXED  AND  LAMINAR  FLAME 


W 


1.^.1 ,  i  '.Ul .  .1  .l,m  I;  J  j^-  ,-^t  l^fllJiU1 
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t  =  Ay  e*KES/RT^)  “‘  111-1 

The  notations  used  here  are  explained  in  Chapter  XL 

With  the  above  assumptions,  the  energy  balance  at  the  surface  is 
given  by  : 

Aj  =  Sp  r  [«,  CT,  -T„)  +w  qap  +0-W)  q4] 

...  m-2 

If  the  thickness  of  the  reaction  zone,  L,  is  specified  as  a  function 
of  pressure  (P),  flame  temperature  (T_),  mixture  ratio  (w),  etc. , 
then  Equations  ni-1  and  III-2  can  be  solved  for  burping  rates  as 
functions  of  P  and  w. 

The  thickness  of  the  reaction  zone,  L,  has  been  estimated  as 
follows  : 

For  a  second-order  reaction,  the  reaction  rate  can  be  written  as 


reaction  fates  B  ?g  w(l-W)  Sx J,(- E*/*Tf) 


...  ra-3 

From  mass  continuity,  we  get 

r  !p--  L  B  C  w(i-w)  e«(.(-EywTf) 

...  Ill -4 


-v' 
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It  should  be  noted  here  that  the  reaction  zone  thickness,  L, 

considered  is  similar  to  L„  considered  in  the  GDF  model, 

2,  ch 

except  that  a  term,  w  (1-w),  is  included  in  the  present  model  to 
account  for  the  dependence  of  reaction  rate  on  the  concentrations 
of  fuel  and  oxidant 

Substituting  the  expression  for  L,  given  by  Equation  III-4  in 
Equation  XII-2,  the  following  equation  for  the  burning  rate,  r,  has 
been  obtained  : 


Equation  III-l  and  ITI-2  were  solved  numerically  for  r  as  a 
function  cf  w  at  various  pressures.  The  input  values  used  are 
listed  in  Table  III-l.  The  results  are  plotted  in  Figure  III-3  and 
are  discussed  in  the  following  section. 

III.  B.  4  DiscussJ-on 

In  Figure  IIT-3,  the  burning  rates,  r,  predicted  by  the  simplified 
model  for  a  polyester/ AP  propellant  are  plotted  against  mixture 
ratio,  w,  for  different  pressures  (P),  and  these  r  vs.  w  curves  are 
compared  with  a  plot  of  flame  temperature,  T„  vs.  w. 


Ag 

= 

2.  OxlO”4 

cal/°C-sec-cm 

cg 

= 

0.3 

cal/°C-gm 

?AP 

1.95 

gm/cc 

9t 

- 

1.27 

gm/cc 

^Ap 

480 

cal/gm 

Qf 

= 

175 

cal/gm 

A 

780 

cm/sec 

s 

E 

s 

15 

Kcal/mole 

e2 

- 

20 

Kcal/mole 

T 

~ 

300 

0 

K 

u 

c t 

-6 

3.  75x10 

* 

T2 

2000°K 

(w=  0.  70) 

T2 

- 

2450°K 

(w:0.  75) 

T2 

- 

2800°K 

(wsrO.  80) 

T~ 

z 

3030°K 

(xu~  0  RfU 
\  —  *  *  *  / 

T 

l2 

mm 

2960°K 

(w=0.  90) 

T 

2 

- 

2500°K 

(wp-O.  S5) 

Calculated  values  assuming  equilibrium  among 
combustion  products _ _ 


TABLE  III-l  LISTING  OF  INPUT  VALUES  USED  FOR 


PREMIXED -LAMINAR  FLAME  MODEL 


BURNING  RATS  (CM/SSC)  Tp  FLAME  TEMPERATURE(°K) 
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The  following  observations  are  made  : 

(i)  The  predicted  burning  rates  attain  a  maximum  value  at  a 

mixture  ratio ,  wr^maxj>  which  is  less  than  the  stoichiometrically 

correct  mixture  ratio,  w  .. 

’  st 

(ii)  The  pressure  does  not  seem  to  affect  the  positions  of  the 

mixture  ratio,  ^  's,  corresoonding  to  maximum  burning  rates. 

(iii)  The  flame  temperature,  T^,,  however,  has  a  maximum 

r 

value  at  a  mixture  ratio  near  the  stoichiometrically  correct  mixture 
ratio,  w  ,  independent  of  pressure. 

fa  t 

These  observations  imply  that  although  the  premixed  flame  model 

can  predict  the  maximum  burning  rates  at  mixture  ratios  different 

from  w  .,  it  is  unable  to  predict  the  shifts  in  w  .  „  with  pressure 

st  r(max) 

variations  such  as  are  observed  in  experimental  investigations. 
Hence,  two  possibilities  were  considered. 

(i)  Due  to  nonequilibrium  chemical  effects,  actual  flame 
temneratures  mav  be  different  from  fhr>c;p  rnlrulntpd  hv  pwiminor 


chemical  equilibrium  among  the  combustion  products,  or 

(ii)  The  surface  of  the  solid  propellant  may  not  be  driven  by  the 
heat  conducted  from  a  single  premixed,  laminar  flame,  as  assumed 
in  the  simplified  model,  but  by  multiple  flames  occuring  in  the  gas 
phase. 
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III.  C.  FLAME  TEMPERATURE  MEASUREMENT 


In  order  to  check  the  possibility  of  non  equilibrium  chemical 
effects,  flame  temperatures  of  loose-granule  burners  were 
measured  at  atmospheric  conditions,  for  various  mixture  ratios. 
The  measured  flame  temperatures  were  then  compared  with 
theoretical  flame  temperatures  calculated  by  assuming  equili¬ 
brium  among  the  combustion  products. 

A  modified  Sodium-D  line  (Na-D  line)  reversal  technique 
was  found  suitable  for  flame  temperature  measurements.  Since 
the  technique  requires  the  flame  to  be  optically  accessible, 
loose-granule  burners  were  fabricated  by  compacting  mixtures 
of  AP  and  polystyrene  in  lead  foil  tubes  rather  than  into  stain¬ 
less  tubes  or  glass  tubes  as  had  been  used  previously  (10). 


T'L  a  f  Ui  n 

ino 


of  the  lead  foil  time  was  chosen  sn  that  the  re¬ 


gressing  surface  of  the  mixture  and  the  melting  rate  of  the  lead 
foil  were  the  same.  The  possibility  of  chemical  effects  of  lead 
on  the  flame  was  checked  by  measuring  the  burning  rate,  which 
is  very  sensitive  to  the  flame  temperature,  for  loose-granule 
burners  in  lead  and  in  stainless  steel  tubes.  These  two  burning 


Gl 

rates  agreed  within  two  percent.  The  thermal  effects  of  the  lead 
tube  were  also  determined  to  be  negligible  since  the  heat  required 
to  melt  the  lead  is  considerably  less  than  that  required  to  pyrolyse 
AP  and  fuel  (e.  g,  480  cai/gm  for  AP,  175  cal/grn  for  a  typical 
fuel,  and  G  cai/gm  for  lead). 

The  details  of  the  modilied  Na-D  line  reversal  technique,  descrip¬ 
tion  of  the  apparatus  used,  and  the  errors  involved  in  the  raeasui  u- 
ments  are  discussed  in  Appendix  A. 

For  comparison  with  experimental  values,  the  adiabatic  flame 
temperatures  of  AP/polystyrene  systems  were  calculated  for  various 
mixture  ratios  using  a  computer  program  developed  by  NASA  (22) 
and  modified  for  use  on  Stevens'  PDP  10  computer. 

III.  C.  1  Results  and  Discus  si  o_n 

Experimentally  measured  flame  temperatures  are  compared 

mini  r.nllxr  /in  Inn  IntnrJ  H  o  m  r»  fomnnv'nfnroc  in  TTi  oiir^  TTT_4 

f  *  ll>It  ilMi/i  V/l  i.  vuu  J  VkV4  V>VtAUVV/U  4.A  ^  vs/  «*»  *  v  S4  s  *♦ 

The  maximum  flame  temperature  in  both  cases  occurs  at  the  same 
mixture  ratio.  This  excludes  the  possibility  of  nonequilibrium 
chemical  effects  being  responsible  for  the  maximum  burning  rate 
occuring  at  mixture  ratios  away  from  w 

Hence,  the  second  possibility,  that  of  multiple  flames  in  the 
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FIGURE  IIM  ; 
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gas  phase  which  supply  the  heat  energy  needed  to  pyrolyse  the  solid 
surface  of  the  propellant,  has  been  considered  to  be  a  reasonable 
one. 

As  a  first  step,  the  comprehensive  models  which  include  more 
than  one  heat  source  were  investigated  further. 
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CHAPTER  IV  -  NUMERICAL  DIAGNOSTIC  TESTS  MADE 

ON  COMPREHENSIVE  MODELS 

IV.  A  INTRODUCTION 

In  the  previous  Chapters  (II  and  in),  the  reasons  were  stated  for 
making  numerical  tests  on  the  existing  comprehensive  models  of  AP 
composite  propellant  burning.  These  many  such  tests  were  not  made 
by  the  originators  of  these  models.  Therfore,  as  part  of  the  present 
investigation,  a  number  of  numerical  tests  were  made  on  the  GDFC 
model,  the  GDFD  model,  the  HR  model,  and  the  MF  model. 

IV.  B  GDFC  MODEL  TESTS 

The  following  tests  have  been  made  on  the  GDFC  model  : 

(i)  A  comparison  of  the  GDFC-model  results  with  experimen¬ 
tal  data  (r  vs.  P). 

(ii)  An  investigation  of  the  influences  of  different  matching 
points  on  the  predicted  results  of  the  GDFC  model. 

(iii)  An  investigation  of  the  influence  of  mixture  ratio  on  the 
burning  rate  as  predicted  by  the  GDFC  model. 
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IV.  B.  1  Comparison  with  Experimental  Bata 

While  presenting  predicted  results  of  the  GDFC  model,  Steinz 

(2)  did  not  compare  them  witn  any  experimental  data.  Hence,  for  the 

purpose  of  comparison,  a  propellant  (AP/LP3,  dQXs78}*.,  and 

w  s  0. 65)  for  which  data  have  been  found  to  correlate  well  with 

Sumrnerfield  equation  (2,  6)  isee  Figure  n-7)  was  chosen*  Following 

Steinz  (2),  a  matching  point  at  Pc  2, 04  atm.  was  chosen  to  determine 

the  emoirical  constants  A  ,  o(,  and  fc.  Using  the  method  of 

s  ‘ 

solutions  described  in  Figure  11-12.  GDFC  equations  were  solved 
for  the  burning  rate.  r.  as  a  function  of  P.  The  input  values  used 
are  given  in  Table  IV-1. 

The  results  of  the  GDFC  model,  using  Steinz 's  matching  point 
("iow-pressure  matching";  .  are  compared  in  Figure  IV-1.  It  is 
found  that  the  GDFC  model  predicts  lower  burning  rates  above 
2,  04  atm.  and  higher  rates  below  2.  04  atm.  The  error  in  correla¬ 
tion  between  the  model  and  experimental  data  is  about  30%  at  luO  atm. 

IV.  B.  2  J^encejrf  ^5atching__Points 

In  order  to  study  the  influence  of  matching  points  on  the  predic¬ 
tions  of  the  GDFC  model,  two  matching  points  other  than  that  used 
by  Steinz  (13.6  atm.  and  100  atm.  )  were  tested.  The  empirical 

constants  A  ,  o<  .  and  ft.  were  determined  at  these  matching  points 
s 


2.0x10  cal/°C-cm~sec 

0.3  cal/°C-cm-sec 

1. 95  gm/'cc 

1. 27  gm/'cc 

480  cal/  gm 

175  cal/gm 

810  cal/gm 

300  °K 

78  microns 

0.65 

2351  °K 

15  Kcal/mole 

20  Kcal/  mole 


MATCHING  POINT 

P  (ATM) 

(°X) 

,As 

(cm/ sec) 

- -  ■  L  ■ 

<X 

2,04 

900 

780 

1. 161xl0"5 

13.6 

1000 

920 

1.067X10'3 

100.0 

1100 

1073 

8.  32xl0“6 

TABLE  IV- 1  INPUT  VALUES  USED  FOR  "GDFC”  AND  ’  GDFD 
MODELS 
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(see  Table  IV-ly,  and  the  GDFC  equations  were  numerically  solved. 
These  results  also  are  plotted  in  Figure  IV-1,  which  show  that  for 
the  high-pressure  matching  point  (100  atm. )  the  predicted  burning 
rates  are  higher  than  tne  exaerimentai  values,  and  the  perr  entage 
error  in  r  at  2„ 04  atm.  is  about  dQV'g.  For  a  mid-pressur  matching 
point  (13.6  atm. ).  the  GDFC  results  were  in  considerable  disagree¬ 
ment  with  experimental  data  „ 

IV.  B.  3  InfLu ence  of  Mixtur e_ Ratio 

A  further  test  involving  the  mixture  ratio  influence  on  r  as 

predicted  by  the  GDFC  model  was  considered.  Numerical  solutions 

to  the  GDFC  equations  were  obtained  by  considering  a  typical 

propellant  whicn  was  considered  earlier  by  Steinz  (2).  Steinz’s 

values  for  the  constants  A  .  and  A  were  used  for  calculations 

1  s 

at  otner  mixture  ratios  (w  ranging  from  0. 7  to  0,  95),  The  flame 
temperatures  used  as  input  values  are  listed  in  Table  ffl-1.  The 
resulting  predictions  of  the  GDFC  model  and  expei  imental  data  are 
compared  ir.  Figure  IV-2.  The  following  conclusions  can  be  drawn  : 

(i)  r  increases  monotonicaliy  witn  w.  ana  no  maximum  burning 
rate  is  predicted  by  the  GDFC  model,  This  is  a  consequence 
of  the  collapsed  AP  flame  assumption,  which  allows  more 
heat  to  oe  generated  at  the  surface  as  w  is  increased. 
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(iil  The  mixture-ratio  sensitivity,  djsT,  predicted  by  the  model 

'  iw 

is  of  the  same  order  of  magnitude  as  those  shown  by 
experimental  data  at  low  pressures,  but  at  higher  pressures, 
it  is  too  low. 

From  these  tests,  it  was  concluded  that  the  GDFC  model  is  not  a 
very  reasonable  model.  Therefore,  it  was  not  considered  for  any 
further  tests  or  modifications. 

IV.  C  GDFD  MODEL  TESTS 

The  GDFD  model  was  also  subjected  to  the  same  tests  which 
were  made  on  the  GDFC  model.  Figure  1^-3  shows  a  comparison 
between  GDFD  results  and  expel  im-ntal  data  for  various  matching 
points,.  The  matching  points  used  in  the  GDFD-model  calculations 
were  obtained  from  the  GDFC-model  calculations.  Since  the  GDFD- 
model  predictions  of  r  are  always  less  than  the  corresponding  GDFC- 
model  predictions,  burning  rates  given  by  the  experimental  data  and 
those  predicted  by  the  GDFD  model  are  not  the  same  at  various 
matching  points.  It  is  clear  from  Figure  IV-3  that  the  disagreement 
between  GDFD  results  and  experimenta!  data  is  great  at  low 
pressures,  whe'e  the  distended  AP  flame  assumption  is  supposed  to 
be  reasonable.  Hence,  any  modification  to  be  made  ori  the  GDFD 
model  should  be  aimed  at  improving  it  in  the  low  pressure  range. 


jfiiiiltt  iTtf ‘rtiraii  ii  r  i  f  iiir  i  nr 
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IV.C.l 


Influence  of  Mixture  Ratio 


The  results  of  the  GDFD  model  with  regard  to  the  mixture-ratio 


influences  on  the  burning  rate  are  shown  in  Figure  IV-4.  T^e 


following  observations  are  made  : 


(i)  r  vs.  w  results  predicted  by  the  GDFD  model  attain  a 


maximum  at  w's  different  from  w  t  at  low  pressures.  This 


trend  is  in  accordance  with  those  shown  by  experiments. 


(ii)  The  mixture  ratio-sensitivities,  predicted  by  the  GDFD 


model  are  much  lower  than  those  observed  experimentally. 


Because  of  the  observed  qualitative  agreement  between  predicted 


and  experimental  trends,  a  few  more  numerical  diagnostic  tests 


were  made.  The  purpose  of  these  tests  was  to  establish  a  basis  for 


modifying  the  model  in  order  to  increase  the  too-low  mixture-ratio 


sensitivity  of  the  model,  while  still  maintaining  its  qualitative 


correct  predicted  trends. 


IV.C.l. a  Modified  "a”  and  "b” 


Two  possible  sources  of  mixture-ratio  effects,  which  do  not 


appear  in  the  GDFD  model  were  introduced  by  modifying  the  GDFD 


model  expressions  for  "a"  and  ,fb"  (Equations  11-13  and  XI-14). 


Though  not  accounted  for  in  the  original  GDF  model  formulation 


(2,.  tne  chemical-time  parameter  "a"  ,  is  sensitive  to  mixture-ratio 
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variations  because  "a"  depends  on  the  reaction  rate  of  the  granular 
diffusion  flame.  In  the  original  formulation  of  the  GDFD  model , 
the  reaction  rate  was  expressed  in  terms  of  the  flame  temperature, 
pressure,  and  other  parameters  (see  Equation  II-6),  Because  the 
reaction  rate  is  expected  to  depend  on  the  concentration  of  the 
reactants,  a  modified  expression  for  the  reaction  rate,  given  by 
Equation  HI-4,  has  been  considered  in  the  present  investigation. 
Denoting  the  modified  a  by  a',  the  following  relation  can  be 
established  between  a  and  a’  : 


CL  *** 


...  IV-1 


Since  the  empirical  constant,^  ,  varies  directly  with  a,  a  modified 

i 

©(  ,  viz. ,  c(  can  be  written  as  : 


o( 


•  •  • 


IV-2 


The  diffusion-time  parameter.  ”b",  considered  in  the  original 
GDF  model  was  assumed  to  be  sensitive  to  dQX,  only,  ance  jA.  , 
the  mass  of  each  fuel  pocket,  depended  only  on  dQX.  It  is  reasona¬ 
ble  to  assume  that  the  mass  of  each  fuel  pocket  can  vary  with  the 
mixture  ratio  also,  is  expected  to  decrease  as  w  increases.  The 
following  functional  relationship  has  been  considered  between  Jw. 
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and  w  : 


|a  ~  (i-w)  f(dox) 
*s  & 

Since,  b^p  ,  biw(i-w£  as  b~f$  ,  the 

ft  ~  ft  0»w) 


then 

V* 


...  rv-s 


...  IV-4 


The  GDFD  equations  (Equations  II-4,  11-17,  and  11-18)  were 
solved  with  both  old  and  new  empirical  constants  viz. ,«(  ,  ,  4 

and^1.  The  following  values  were  used  : 

e<-  g.£x  iok  ;  ft  «  O’  10  99  (Calculated  from  data  of  Ref.  1) 

1  •»it  i 

©1  =  5-T’XIO  ]  ft  -  0* (From  Equations  IV-2  and  IV-4) 

The  predicted  results  of  the  GDFD  model  for  both  new  and  old 
parameters  <<  and  ft  are  compared  in  Figure  IV-5.  Although  the 
trends  of  the  r  vs.  w  curves  have  been  preserved,  the  modification 
has  not  improved  the  mixture-ratio  sensitivity,  ^DJ*. 

IV.  C.  1.  b  Modified  "Pyrolysis  Rate  Equation" 


In  the  original  GDFD  model,  an  overall  surface-pyrolysis  rate 
equation  was  assumed  to  describe  the  surface  pyrolysis  (Equation  n- 
4).  This  equation  does  not  contain  any  term  involving  the  mixture 
ratio,  w.  In  the  MF  model  (19,  20),  a  different  surface  pyrolysis 


ORIGINAL  ’a*  AND  ’V 
MODIFIED  'a*  AND  ’b* 


P  -  1  ATM 
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rate  equation  containing  mixture  ratio,  w,  was  considered 
(Equation  11-22),  Hence,  Equation  II-4  was  replaced  by  Equation  II- 
22  and  the  resulting  equations  of  the  GDFD  model  (Equations  H-17 
and  n-18)  were  solved  for  the  burning  rate  while  varying  the  mixture 
ratio. 

The  results  obtained  using  both  the  original  and  modified  pyroly¬ 
sis  -  rate  laws  are  shown  in  Figure  IV-6.  This  modification  basica¬ 
lly  changes  only  the  magnitude  of  burning  rates  while  preserving 
the  trends  of  r  vs.  w  curves.  It  has  not  increased  the  predicted 
mixture-ratio  sensitivity, 

dw 

IV.  C .  2  Further_  Numerical  _Dj.agnostic  _Tests  (Inver s_e 
Procedure) 


As  the  previous  tests  failed  to  improve  the  GDFD  model,  with 
respect  to  mixture-ratio  sensitivity,  the  following  diagnostic  tests 


WUX't:  made  to  isulate  particular  ill 6  C  1’nt  iii  5 ill a  WuiCu,  when  included, 


might  improve  the  GDFD  model.  It  was  observed  that  both  the 
energetics  (heat  generated  at  the  first-stage  flame,  and  kinetics 
(reaction  time  of  the  first-stage  flame,  TTj )  play  important  roles  in 
the  predictions  of  the  GDFD  model.  These  roles  were  investigated 
by  means  of  an  "Inverse  Procedure". 


ir ’Hiffiffiiiiirn  • 
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The  inverse  procedure  involves  using  experimental  burning-rate 
data  for  various  pressures,  particle  sizes,  and  mixture  ratios  as 
input  values  in  the  GDFD  equations  to  calculate  the  heat  generated  at 
the  first-stage  flame,  AHp  and  the  reaction  time  of  the  first-stage 
flame  tj. 

The  experimental  data  used  here  were  taken  from  Reference  1 
for  AP/ P13  propellants  with  w«0. 7,  and  dQXs  16  }*  and  120 
The  GDFD  model  was  forced  to  predict  r  vs.  w  values  for  ws  0. 75 
and  d^  16  ju  and  120  /<-  by  suitable  choices  of  o^and  .  Such  a 
method  had  to  be  used  since  the  GD FD-mod el-pr edicted-values  were 
lower  than  the  experimental  data  (see  Figure  IV- 3).  This  selective 
choosing  of  values  for  *  and  can  be  justified:  The  purpose  of  the 
inverse  procedure  is  not  to  obtain  quantitative  results  but  rather 
qualitative  results  which  can  provide  clues  for  improving  the  model. 

The  set  of  values  used  for<*  and/3  so  obtained  for  w  *  0. 75  and 


dQ=r  and  120  were  used  along  wim  x  Vo.  it*  iOi  w  —v.o 
in  the  GDFD  equations  to  calculate  and  . 

In  calculating  ,  T(was  assumed  to  be  a  constant  with  respect 
to  d  and  w  but  was  varied  inversely  with  P,  as  given  by  Equation 

UA 

n-16.  The  results  of  the  inverse  procedure  for  calculating  <JIij  are 
shown  in  Figure  IV-7.  These  results  indicate  that  if  the  heat  genera¬ 
ted  at  the  first-stage  flame  were  higher,  depending  onP,  d  ,  and 

UA 
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w,  the  GDFD  model  would  have  predicted  reasonable  .results. 

The  inverse  procedure  was  also  used  to  calculate  the  reaction 
time  of  the  first- stage  flame,  X t ,  from  the  GDFD  equations  by 
assuming4H^  to  be  constant  with  respect  to  P  but  by  assuming  a 
linear  variation  with  respect  to  w  (as  assumed  in  the  GDFD  model). 
These  calculations  yielded  unrealistic  values  for  X ,  (sometimes 
neagative  values  were  obtained).  Therefore,  it  was  decided  to  vary 
the  reaction  time  for  the  first-stage  flame,  X. ,  in  a  reasonable  way 
with  respect  to  w  and  to  solve  the  GDFD  equations  with  the  modified 

.  It  was  noted  that  from  the  results  reported  in  Reference  2 
that  the  reaction  time  for  the  second-stage  flame  in  the  GDFD  model 
viz. ,  X^,  varied  with  both  w  and  P.  Tis  given  by 


PM  C£il 

l  p 


Pvj  T% 
$ 


i 


...  IV-5 


where  M  is  the  average  molecular  weight  of  the  gases,  and  R 
is  the  universal  gas  constant.  All  the  other  quantities  in  Equation 
IV-5  have  the  same  meaning  as  those  described  in  Chapter  II. 

At  low  pressures,  it  was  aibu  observed  that  T^y/as  about  10  times 
greater  than  %  .  Hence,  in  order  to  vary  'tjwith  both  P  and  w,  the 
following  relation  between  t]and  T2was  considered. 
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The  GDFD  model  equations  were  solved  numerically  with  %  - 
and  The  results  are  plotted  as  r  vs.  w  for  P  1  atm.  and 

are  compared  in  figure  IV-8.  It  is  clearly  seen  that  the  r  vs.  w 
curves  corresponding  to  tjso-il^have  a  steeper  slope  than  the  r  vs.  w 
corresponding  to  These  results  show  that  if  the  kinetics  of 

the  first-stage  flame  can  also  be  varied  by  some  mechanism,  then 
the  GDFD  model  can  predict  increased  mixture-ratio  sensitivities, 

iw 

The  results  of  the  inverse  procedure  clearly  indicate  that  the 
GDFD  model  can  be  improved  if  the  first-stage  flame  were  perturbed 
by  some  mechanism  vhich  would  change  the  heat  generated  at  the 
first-stage  flame  and  the  reaction  time  of  the  first-stage  flame. 

IV.  D  HR  MODEL  TESTS 

The  heterogeneous  reaction  model  (HR)  has  been  tested  in  a 
comprehensive  manner  previously  (17. 18).  Thesp  tests  showed 
the  model  to  prelict  burning  rate  behaviour,  which  is  qualitatively 
in  agreement  with  experimental  results.  However,  in  these  tests, 
mixture  ratio  influences  were  considered  only  for  a  narrow  range 
of  variations  in  w  (0.6  to  0.7).  Hence,  the  trends  of  r  vs.  w  curves 
were  not  evident  in  these  numerical  model  tests.  Therefore,  as 
part  of  the  present  investigation,  a  further  test  involving  mixture- 
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ratio  influences  was  made, 

I V.  D .  1  Influence _of  _Mixture  _Rjitio_ 

The  equations  of  the  HR  model  were  numerically  solved  for 
mixture  ratios  beyond  wsO„  7.  The  program  for  numerical  calcula¬ 
tions  was  supplied  by  the  model's  originator.  The  basic  equations 
employed  in  the  HR  model,  viz. ,  mass  conservation  and  energy 
conservation,  allow  the  model  to  be  extended  to  the  limiting  case 
when  w:  1.  0  (pure  AP).  However,  an  error  developed  in  the 
program  prevented  the  calculations  from  being  carried  beyond 

ws  0.95  for  d  =  20k  and  wsO,  9  fo"  d  *  200J*.  Hence,  the 
ox  *  ox  f  ’ 

numerical  calculations  were  not  carried  beyond  ws  0. 95  for  d  ox 

and  ws  0.9  for  d  =  200k. 

ox  * 

The  results  of  the  HR  model  are  plotted  as  r  vs.  w  in  Figure  rV  9 
for  varying  P  and  dQx.  It  can  be  seen  clearly  that  the  HR  model  does 
not  predict  maximum  values  for  r.  This  trend  is  contradictory  to 
the  experimental  data  trends. 

I V,  D .  2  Pis£uj?  sion 

The  major  objection  to  considering  the  HR  model  any  further 
is  that  the  heterogeneous  reaction  assumed  in  the  model  seems  to  be 
an  important  mechanism  which  influences  the  burning  rate  prediction 


w,  WEIGHT  FRACTIOI!  OF  OXIDANT 

FIGURE  IV-9  INFLUENCE  OF  i  IIXTURE  RATIO  AS  PREDICTED  BY  THE 
HEIEROGENECUS  REACTION  MODEL 
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of  this  model.  However,  no  experimental  evidence  has  been  found  so 
far  to  support  the  existence  of  such  reaction.  This  fact  coupled  with 
the  fact  that  the  model  predicts  no  maximum  for  the  burning  rate,  led 
toward  abandonment  of  this  model  as  a  basis  for  any  modifications  or 
additional  tests. 

IV.  E  MF  MODEL  TESTS 

Numerical  solutions  were  obtained  for  the  multiple  flame  model 
(MF)  equations  by  using  a  computer  program  supplied  by  the  origina¬ 
tors  of  the  MF  model.  The  complete  set  of  input  data  found  in 
Reference  19  were  used,  and  the  results  reported  in  Reference  19 
were  reproduced.  However,  numerical  calculations  made  with  the 
incomplete  set  of  input  data  found  in  Reference  20  and  additional  data 
suggested  by  an  originator  of  the  MF  model  (23)  yielded  results  which 
did  not  agree  with  those  i  epcrted  in  Reference  20.  In  Figure  IV-10 
the  results  obtained  presently  and  those  reported  in  Reference  20  are 
shown.  Two  possibilities  may  be  considered  for  the  discrepency 
shown  : 

(i)  an  error  might  have  developed  in  the  program,  or 
(ii)  a  still  different  set  of  input  data  might  have  been 
used  to  get  the  results  reported  in  Refernce  20. 

In  any  case,  the  results  of  Reference  20  should  not  be  considered 
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FIGURE  Vi-20  :  RESULTS  OF  HE  MULTIPLE  FUME  ITOEL 


useful  until  a  complete  set  of  input  data  is  available. 

Besides  the  above  mentioned  discrepency,  there  are  a  few  other 
demerits  of  the  MF  model. 

In  estimating  the  diffusion  flame  heights,  the  analysis  of  Burke 
and  Schumann  was  used  (22).  Since  the  flame  structure  assumed  in 
this  model  requires  all  the  flames  to  envelope  the  oxidant  particle, 
only  fuel-rich  diffusion  flames  ("underventilated")  can  be  considered. 
This  restricts  the  propellant  composition  that  can  be  described  by  the 
MF  model  to  fuel-rich  compositions.  Hence,  mixture-ratio 
variations  cannot  be  considered  over  a  wide  range. 

Another  major  drawback  of  this  model  is  the  assumption  that  the 
primai  y  flame  temperature  is  equal  to  the  adiabatic  flame  tempera¬ 
ture  of  the  propellant.  The  primary  flame  is  the  reaction  between  a 
fraction  of  AP  products  and  a  fraction  of  fuel  vapors.  The  fraction 
of  AP  products  that  react  in  the  primary  flame,  A  >  is  a  function 

i  p 

of  pressure  and  particle  size.  Therefore,  the  primary  flame  tempe¬ 
rature  should  depend  on  these  variables.  The  adiabatic  flame  tempe¬ 
rature  does  not,  however,  depend  on  either  the  pressure  or  the 
oxidant  particle  size.  Thus  the  validity  of  assuming  the  calculated 
adiabatic  flame  temperature  and  primary  flame  temperature  as 
equal  is  questionable. 


WSZ£  -w.  ** '  tjh&zgm 


For  the  reasons  stated  above,  the  MF  model  was  not  considered 
for  any  further  modifications  or  tests. 

IV. F  CONCLUSIONS 

The  results  of  numerical  diagnostic  tests  made  on  several 
comprehensive  models  have  indicated  the  GDFD  model  to  be  a 
promising  one,  i.  e. ,  one  which  might  possibily  be  modified  to 
predict  results  in  better  accord  with  experimental  data.  Such 
improvement  should  be  in  two  categories  : 

(i)  an  increase  in  the  predicted  burning  rates  at  lower 
pressures,  and 

(ii)  an  increase  in  the  mixture-ratio  sensitivity  of  the 
burning  rates. 

The  results  obtained  by  the  ‘'inverse  procedure"  indicated  that 
the  GDFD  model  could  be  substantially  improved  if  the  first-stage 
flame  were  perturbed.  A  mechanism  for  such  perturbation  is  that 
of  intermixing  of  fuel  vapors  into  the  AP  decomposition  flame  and 
of  consequent  additional  reaction.  This  mechanism  was  long  ago 
suggested  by  Chaiken  and  Anderson  (19, 20)  and  quite  recently 
accounted  for  in  the  MF  model  (19,  20). 

A  modified  granular  diffusion  flame  model  incorporating  the 
above  mentioned  has  been  formulated  and  tested  in  the  present  study. 
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CHAPTER  V  -  MODIFIED  GRANULAR  DIFFUSION  FLAME  MODEL 


A  modified  granular  diffusion  flame  model  (MGDF),  incorporating 
the  intermixing  of  fuel-vapors  into  the  AP  decomposition  flame  and 
the  consequent  additional  reaction  with  the  AP  decomposition  flame, 
is  described  below.  The  major  assumptions  made,  the  resulting 
equations,  and  the  numerical  tests  made  on  the  MGDF  model  to  assess 
its  validity  are  also  described  in  the  following  sections. 


V.  A  DESCRIPTION  OF  THE  MGDF  MODEL 


The  MGDF  model  considers  a  two-stage  flame  in  the  gas  phase  to 
provide  the  necessary  heat  required  to  gasify  the  surface.  The  first- 
stage  flame  is  a  perturbed  AP  decomposition  flame  and  the  second - 
stage  flame  is  a  granular  diffusion  flame.  The  perturbation  of  the 
AP  flame  is  due  to  the  intermixing  of  fuel  vapors  into  the  AP  flame 


and  the  additional  reaction  with  the  Af  flame. 


jfigure  V-I  describes 


the  MGDF  model. 


V.  B  MAJOR  ASSUMPTIONS 

The  following  assumptions  are  made  : 
(i)  One  dimensional. 

(ii)  Steady  state. 
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(iii)  Specific  physico-chemical  processes  : 

The  physical  processes  considered  are  heat  conduction 
and  granular  diffusion. 

The  chemical  reactions  considered  are  fuel  pyrolysis, 
dissociative  sublimation  of  AP  into  ammonia  and 
perchloric  acid,  a  perturbed  AP  decomposition  flame, 
and  a  granular  diffusion  flame. 

(iv)  Linear  temperature  distribution  in  the  gas  phase. 

(v)  Uniform  thermal  properties. 

V.  C  EQUATIONS 

As  a  consequence  of  the  above  mentioned  assumptions  and 
modification,  the  resulting  equations  for  the  MGDF  model  are  the 
same  as  those  for  the  GDFD  model  (Equations  II-4,  11-17,  and  11-18) 
Only  the  input  values  related  to  the  first-stage  flame,  viz. ,  AHy 
the  heat  generated  at  the  first-stage  flame,  and  'Cj ,  the  reaction  time 
of  the  first-stage  flame,  are  different.  In  the  GDFD  model  these 
two  input  values  are  determined  from  the  deflagration  data  for 
pure  AP,  but  in  the  MGDF  model,  they  are  determined  from  the 
intermixing  mechanism  of  fuel  vapors  as  functions  of  pressure, 
oxidant  particle  size,  and  weight  fraction  of  oxidant,  as  described 


below. 
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Following  the  notations  used  in  Chapter  II,  the  set  of  equations 
to  be  solved  in  the  MGDF  model  are  : 

r-  as  e*fc>(-E5/frrs) 


...  v-i 


T  -- 


...  V-2 


r  =  U  (■ 


A,(Tt-T„) 


1 


Cgrn-To)  +w  QaP  +U-W)  c^-w  AH,  j 


...  V-3 

The  modified  input  values  and  tjwere  determined  as  follows. 
* 

Let  T  denote  the  true  temperature  of  the  perturbed  AP  flame. 

Due  to  the  additional  reaction  between  the  fuel  vapors  that  intermixed 

* 

with  the  AP  flame,  could  be  determined  if  the  amount  of  fuel 

vapor  that  reacted  with  the  AP  products  could  be  specified.  Once 
$ 

has  been  determined,  then  following  the  original  GDFD  model, 
the  following  equation  could  be  used  to  determine  ; 


AH,  =  S<Wo)+  <=* 


V-4 


A? 


Assuming  that  the  perturbed  AP  flame  is  a  second  order  reaction, 
is  given  by; 


or 


A/ 


I 


?3  -E/Tuf) 


•  <  « 


V-5 


The  constant  Ci  is  determined  from  the  condition  that,  for 

%  ^  (f 

Tl"TAP'  (Equation  11-16). 

-12  sec -atm 

C1  was  found  to  be  1.  8  x  10  ( — - - ). 

1  °K 

In  order  to  determine  the  amount  of  fuel  vapor  that  intermixes 

* 

with  the  AP  flame  and  hence,  ,  the  following  model  was 
considered.  Following  Figure  V-l,  let  dQx  denote  the  oxidant 
particle  diameter,  Lj  be  the  height  of  the  first-stage  flame  above 
the  surface  of  the  propellant,  if  rh^-  represents  the  mass  flow  rate 
(gm/sec)  of  the  fuel  that  intermixes  through  the  cylindrical  surface 
area  surrounding  the  AP  particles,  then 


=  ? y  V  A 
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where  £  is  the  average  gas  phase  density,  V  is  die  diffusion 
velocity,  and  A  is  che  area  available  for  mixing. 

V  can  be  written  as  : 

V  cr 


2> 

d 


. .  V-7 


os 


.1.75 


where  D  (the  diffusion  coefficient).-^ 


...  V-8 


A  =  Tf  dw  L, 


...  V-9 


Using  Equations  V-7  and  V-9  in  Equation  V-6,  we  get 

riif*  f  T>  l,  •••  v-10 


In  Equation  V-10,  P  and  D  are  related  to  P  and  T  by  the  perfect 


tew  and  by  Equation  V-8,  respectively,  and  is  given  by, 


L'  (-r) 


...  v-ll 


where  m"  i  the  mass  flux  of  the  propellant  (a  $  t* ). 


T  *•'"  v  *P 

The  reaction  time  appearing  in  Equation  V-ll  is  given  by  Equation 


V-5.  Since  is  an  unknown  quantity,  introducing  it  into  the  expre¬ 


ssion  for*lj,  would  involve  an  iterative  method  to  determine  . 
Hence,  for  mathematical  simplicity,  it  was  assumed  that  an  average 


temperature  T  ,  given  by  the  following  equation  be  used  in  the 
&v 


V-12 


place  of  in  Equation  V-ll, 


T  = 

av 


T,+t; 


AP 


2- 


where  T  and  TAO  are  the  propellant  and  AP  flame  temperatures, 
2  Air 

respectively. 

The  reaction  time,  used  in  the  present  model  is  given  by; 

I 

**  ...  V-13 

Equation  V-10  can  be  reduced  to  the  following  form  by  using 
Equations  V-ll  ,  11-12,  and  11-13  : 


Tftv 


a- 'Pi 


ft  -  m"  — - — — 

4  T  V  ^t) 


...  V-14 


Letting  m"  represent  the  mass  flux  of  AP,  then  the  mass 
ox 


flow  rate  m  turn /seel  for  a  simile  narticle  is  given  by  ; 

ox  •*  •  ~ 


•  *  M  .  | 

in  -  bv  rr  ( Sex) 

'ox  Or.  \  2.  / 


...  V-15 


m^ 

The  ratio  ~'m“x  can  then  be  written  as  : 


.  /  .  ** 

J  Mt  -  • 

\rh 


MB 

Tftv 


\  KJ 


Ox  '  ox 


V-16 


The  ratio 


is  given  by  : 


in 


n 

r 

it 


OX 
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Finally  then,  combining  Equation  V-16  and  V-17,  the  ratio  of  fuel  to 
oxidant  in  the  first-stage  flame  is  proportional,  to  functions  of 
pressure,  oxidant  particio  sine,  and  other  propellant  properties  is 
given  by  : 


...  V-18 

C2  is  &  constant  which  can  be  determined  empirically  by  matching 
the  burning  rates  predicted  by  the  model  and  given  by  experiments 
at  a  matching  point.  Details  of  the  determination  of  CL  are  given  in 
Section  V.D. 

Once  the  fuel -to -oxidant  ratio  of  the  first-stage  flame  is  deter - 

*• 

mined  from  Equation  V-18,  can  be  determined  from  thermo- 
chemical  calculations  (22).  It  should  be  noted  that  the  present 
model  predicts  that  the  amount  of  fuel  mixing  with  the  AP  flame  is 
greater  at  low  pressures  and  for  small  oxidant  particle  sizes 
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(Equation  V-18),  which  is  in  agreement  with  the  results  indicated  by 
the  inverse  procedure  in  Chapter  IV  (see  Figure  TV-7). 

V. D  METHOD  OF  SOLUTION 

The  following  steps  were  used  in  solving  the  MGDF-model 
equations  : 

(i)  A  propellant  whose  burning  rate  data  correlated  with 
Summerfield's  equation  was  chosen.  The  empirical 
constants  Ag,  e{,  and  ^  were  determined  as  describ¬ 
ed  in  Chapter  II. 

(ii)  A  pressure  (P)  was  chosen  in  the  low  pressure  region, 

where  the  modification  on  the  GDFD  model  would  be 
mr 

elective.  _  was  determined  from  Equation  V-18 
mox 

by  assuming  a  reasonable  value  for  Cg  and  for  the 

input  values  of  w,  d  .  etc. , 

ox 

(iii)  For  the  mf  calculated  in  step  (ii),  the  correspon- 


calculations. 

(iv)  The  modified  input  values  and  'Cjwere  determined 
from  Equations  V-4  and  V-13,  respectively. 

(v)  With  the  appropriate  input  values  (steps  i  and  ii), 
the  resulting  equations  (Equations  V-3.  to  3)  were 


■  - 


H&iSj 


I: 


4: 


f 


(Vi) 


(vii) 


V.  E 
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solved  numerically  to  obtain  at  the  selected 


pressure  (Pj.  The  listing  of  the  computer  program  used 
is  given  in  Appendix  B. 


The  calculated  r  was  compared  with  r„  if 

MGDF  EXPT 


they  were  different,  a  new  Cg  was  assumed  and  steps  (iii) 
to  (v)  were  repeated  until  agreement  was  reached. 

The  constant  Cg  determined  in  step  (vi)  was  used  to 
calculate  r  as  a  function  of  P,  dQX,  etc. , 


RESULTS  ARB  DISCUSSION 


W 


The  MGDF  model  has  been  subjected  to  the  following  tests  in 
order  to  check  its  predictive  capability  for  predicting  burning  rate 
variations  with  varying  pressure  (P),  oxidant  particle  size  (a), 
mixture  ratio  (w),  and  initial  temperature  (TQ). 

V.  E.  i  _Pr_essureJEf f ecj:s_as _P£G_dict_cd_ ^yJJ  uoJ^GpF  Mod el_ 

The  MGDF-model  equations  were  solved  for  burning  rates  at 
various  pressures  for  a  propellant  whose  r(P)  data  correlated  well 
with  Summcrfield's  equation,  viz.  Equation  II~1.  The  propellant 
chosen  lor  this  calculation  was  an  AP/LP3  propellant  with  w  *0.  65 
andd  *78^t.  The  experimental  data,  i.  e.  r  vs.  P  data,  were 
taken  from  Reference  6. 
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The  necessary  input,  data  needed  for  the  numerical  calculations 
are  listed  in  Table  IV-1,  The  empirical  constants  A  ,  d,  and  (S 
were  obtained  from  a  high  pressure  matching  (P  100  atm. )  proceed- 
ure  described  in  Chapter  II.  The  constant  C^,  introduced  in  the 
MGDF  model,  was  determined  by  the  iterative  procedure  described 


in  the  previous  section  (Section  V.D)  by  matching  the  burning  rate 

predicted  by  the  MGDF  model  with  the  burning  rate  given  by 

experiment,  at  a  low  pressure  (ps  2.  04  atm. ).  The  valu  e  of  C2 

-17 

obtained  by  the  iterative  procedure  was  found  to  be  1.0  xlO 

/  cm2-atm2v 
'  ok2TW* 

Figure  V-2  shows  the  predictions  of  the  MGDF  model  compared 
with  the  results  of  both  the  GDFC  and  the  GDFD  models  and  also  with 
experimental  data. 

The  following  conclusions  can  be  drawn  from  this  comparison  : 

(i)  The  burning  rates  predicted  by  the  MGDF  model  are 
in  good  agreement  with  the  experimental  data  in  the 
pressure  range  of  1  to  100  atm.  The  maximum  error 
in  this  range  between  the  MGDF  results  and  experimen¬ 
tal  results  is  about  10%. 


The  results  of  the  MGDF  model  are  clearly  much 


better  than  those  of  GDFC  and  GDFD  models.  The 
improvement  is  quite  significant  at  the  low  pressure 


'  region,  where  the  GDFC  and  GDFD  models  predict 

higher  (about  60%  above  the  experimental  burning  rate) 
and  lower  (about  40%  below  the  experimental  burning 
rate)burning  rates,  respectively. 

V.  E.  2  Par  tic!  e_  Size_  Effects  as_  Predicted_by  J^^MGDF 

Model_ 

In  order  to  test  the  MGDF  model’s  ability  to  predict  the  effects 
of  oxidant  particle  size  (d  )  on  the  burning  rate,  an  AP/LP3 
propellant  containing  140  j*AP  particle  size  and  65%  of  AP  by  wreight 
was  considered.  The  experimental  data  (r  vs.  P)  for  this  propellant 
have  been  shown  to  correlate  well  with  Summerfi eld’s  equation 
in  References  2  and  6.  It  should  be  noted  that  the  propellant  chosen 
here  has  the  same  composition  as  that  considered  in  the  previous 
section  to  study  pressure  effects,  i.  e. ,  both  are  AP/LP3  propell¬ 
ants  and  w:  0.  65  for  both  of  them.  Only  the  particle  sizes  of  these 
propellants  are  different. 

As  explained  in  Chapter  II,  the  variations  in  the  particle  size 
affects  only  the  diffusion-time  parameter  b,  and  hence,  the 
empirical  constant  Therefore,  for  the  numerical  calculations 
only  has  been  changed  while  keeping  other  constants  ct,  C^,  C2, 
etc. ,  the  same  as  those  cited  for  78f*particle  size  in  the  previous 
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section  (Section  V.  E.  1). 

*  • 

Since  the  constant  C g  used  for  140|*  and  78^  AP  particle  sizes 
were  taken  to  be  the  same,  an  iterative  procedure  was  not  required 
in  this  case.  Thus  the  MGDF  equations  were  directly  solved  for  the 
burning  rate  at  various  pressures. 

In  Figure  V~3,  a  comparison  of  the  MGDF  results  with  both  the 
GDFC  and  the  GDFD  results  and  also  with  experimental  data  is  made. 

As  seen  clearly  the  predicted  results  of  the  MGDF  model  are  in 
good  agreement  with  the  experimental  data  and  are  much  closer  than 
both  the  GDFC  and  the  GDFD  results.  Since  the  propellant  considered 
presently  contains  a  larger  particle  size  than  the  one  considered  in 
studying  the  effects  of  pressure  on  r,  it  is  concluded  that  the  MGDF 
model  not  only  predicts  the  effects  of  pressure  on  the  burning  rate 
correctly  but  also  the  effects  of  oxidant  particle  size  effects  in 
accordance  with  experimental  data. 

V.E.3  Maxtor  e_  RatioJSffecte  as  ^Predicted  Jby_  the 

MGDF  Model_ 

One  of  the  objectives  in  modifying  the  GDFD  model  was  to 

increase  the  too-low  mixture-ratio  sensitivity,  ,  predicted  by 

$vf 

the  GDFD  model.  To  test  the  MGDF  model  in  this  regard,  the  MGDF 
model  equations  were  solved  for  input  values  appropriate  for  an 
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AP/P13  propellant  with  dQX«  16  jut  (1).  The  empirical  constants 
needed  for  the  numerical  calculations,  viz. ,  Ag,  ,  and  C % 

were  determined  for  w=  0. 75  and  these  values  were  ther  used  for 
other  mixture  ratios  (i.  e. ,  w  a  0.  8,  0.  85,  0.  9,  and  0. 95).  These 
values  are  listed  in  Table  V-l. 

The  MGDF-rnodel  results  are  plotted  as  r  vs.  w  in  Figure  V-4 
where  they  are  compared  with  both  the  GDFD-model  results  and  the 
experimental  data.  It  is  seen  that  the  sensitivities  given  by  the  MGDF 
model  are  much  higher  than  those  given  by  the  GDFD  model.  For 
example,  the  mixture-ratio  sensitivity,  4br,  given  by  the  GDFD 
model  is  about  1. 5  and  the  sensitivity  given  by  the  MGDF  model  is 
about  5, 5  for  P-  2. 04  atm.  The  experimental  results  show  an 
average  value  of  about  5. 0  for  mixture  ratio  sensitivity  at  the  same 
pressure. 

Comparison  of  the  MGDF-model  results  with  experimental  data 
could  not  be  made  beyond  v/sO.  8,  because  experimental  data  are  not 
available  at  present.  Hence,  only  a  qualitative  comparison  can  be 
made  between  the  MGDF  model  predictions  and  existing  analog  data 
with  regard  to  the  trends  of  r  vs.  w.  The  MGDF  model  predicts 
maximum  burning  rates  at  different  mixture  ratios  depending  on  the 
pressure.  Such  a  behaviour  of  the  burning  rate  has  been  observed 
in  the  studies  involving  analogs  of  propellants. 
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FLAME  TEMPERATURES 

T0  -  2450  °K  (w=r  0. 75) 

T2  -  2800  °K  (w  ~  0. 80) 

T2  =  3020  °K  (w  =  0.85) 

T2  -  2960  °K  (w  a  0.90) 

T  =  2500  °K  (w  a  0.95) 

2 

EMPIRICAL  CONSTANTS 

X  =  6. 136xl0~5 

P  =  0. 2659 

1 9 

C  =  1.8x10 

_1  -  ,  ..-18 
u.  -  o.  Oxiu 

2 


TABLE  V-l  -  LISTING  OF  INPUT  VALUES  USED  IN  THE 

MGDF-MODEL  CALCULATIONS 


BURNING  RITE  (CM/SBC) 


0.1b  - - 


0.75  0.8  0.85  0.9  0.95 

w .  WEIGHT  FRACTION  OF  OXIDANT 


FIGURE  V-4  :  MIXTURE-RATIO  INFLUENCES  AS  PREDICTED  BY  THE 
MGDF  MODEL 
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V .  E .  4  *  Initial^  Temp  era  tur  e_  Effects  _as  _Pr  edicted_  by_  the 

MGDF  Model_ 

Finally,  the  MGDF  model  lias  been  tested  for  initial-propellant 
temperature  influences  on  the  burning  rate.  The  propellant  consider¬ 
ed  for  numerical  calculations  was  an  AP/LP3  propellant  with 
dox=  I40|iand  wsO.  65.  Initial  temperatures  were  varied  between 
300°K  and  200°K.  The  results  of  the  numerical  calculations  are 
plotted  in  Figure  V-5. 

The  MGDF  model  predicts  that  the  burning  rate  increases  as 
the  initial  temperature  is  increased  over  a  wide  range  of  pressure. 
This  prediction  is  in  qualitative  agreement  with  the  experimental 
data.  The  initial  temperature  sensitivity.^)*^.  given  by  the  MGDF 

K 

model  is  about  0.  08%/oq  ,  whereas  the  experimental  results 
available  for  a  different  propellant  (AP/PBAA)  (7)  is  about 
0.  2%/oq.  The  lower  initial  temperature  sensitivity  given  by  the 
MGDF  model  may  be  due  to  the  fact  different  propellants  were 
considered  in  the  experimental  and  theoretical  investigations.  Hence, 
no  definite  conclusions  can  be  made  about  the  merits  or  demerits  of 
the  MGDF  model  based  on  this  test. 


AP/LP3  PROPELLANT 


INFLUENCES  OF  THE  INITIAL  TEMPERATURE  AS  PREDICTED  BY  DC  MGDF  MODEL 


110 


V.  F  '  CONCLUSIONS 

On  the  basis  of  the  above  tests  made  on  the  MGDF  model,  the 
following  conclusions  are  drawn  : 

(i)  The  MGDF  model  predicts  reasonable  burning  rate 
variations  with  pressure  in  the  range  of  1  to  100  atm. 
The  MGDF  gives  results  which  are  much  closer  to 
experimental  data  than  those  of  the  GDF'*  and  GDFD 
models. 

(ii)  The  oxidant  particle  size  effects  predicted  by  the 
MGDF  mod^l  are  also  in  agreement  with  experimental 
data  and  in  closer  agreement  than  those  of  the  GDFC 
and  GDFD  results. 

(iii)  The  mixture-ratio  sensitivities  predicted  by  the  MGDF 
model  are  essentially  the  same  as  those  given  by 
experimental  data  and  are  much  higher  than  those 
predicted  by  the  GDFD  model. 

(iv)  The  predicted  dependence  of  the  burning  rate  on  the 
initial  propellant  temperature  is  qualitatively  in 
agreement  with  experimental  data. . 


CHAPTER  VI  -  MAJOR  CONCLUSIONS  AND  SUGGESTIONS 


FOR  FUTURE  RESEARCH 

VI.  A  MAJOR  CONCLUSIONS 

The  present  experimental  and  theoretical  investigations  with  a 
loose- granule  burner  and  the  theoretical  investigations  with  a 
premixed,,  laminar  flame  model  indie, *  hat  the  gasification  of  the 

surface  of  a  composite  solid  propellant  is  probably  driven  by  multiple 
heat  sources  (flames)  in  the  gas  phase.  The  influence  of  mixture- 
ratio  on  the  burning  rate  cannot  be  explained  using  simplified 
mechanistic  models  which  consider  only  one  flame  in  the  gas  phase. 

One  of  the  important  reactions  in  the  gas  phase  has  been  shown, 
in  this  study,  to  be  a  reaction  between  fuel  vapors  intermixing  into 
the  AP  flame  and  AP  products  in  the  area  surrounding  each  AP 
particle.  Although  this  mechanism  was  proposed  long  ago  by 
Chaiken  and  Anderson  (14),  its  importance  has  not  been  established 
until  this  investigation. 

A  modified  granular  diffusion  flame  model,  formulated  by 
including  the  above  mentioned  mechanism  in  the  GDFD  model,  gives 
results  which  agree  with  experimental  results  more  closely  than 
do  results  from  the  GDFD  model.  The  MGDF  model  has  been 
comprehensively  tested  against  experimental  data  and  is  claimed 


to  be  the  most  valid  model  in  that  it  is  based  on  reasonable  meehanis- 
* 

tic  grounds  and  it  is  capable  of  predicting  the  influences  of  important 
variables  on  the  burning  rate  as  indicated  by  experimental  results. 

VI.  B  SUGGESTIONS  FOR  FUTURE  RESEARCH 

There  is  a  need  for  further  more  comprehensive  data  for 
composite  solid  propellants,  especially  for  mixture-ratio  (w)  varia¬ 
tions.  Hence,  it  is  suggested  that  future  experimental  investigations 
be  planned  to  obtain  comprehensive  data  at  least  for  those  which  obey 
Summerfield's  equrtion. 

The  simplified  assumption  made  in  the  MGDF  model  with  regard 

to  the  change  in  the  first-stage  flame  kinetics  (*fcp  (i.  e.  Tav  was 

* 

used  in  Equation  V-13  instead  of  T  )  should  be  removed  by  replacing 

* 

Tav  ky  and  an  iterative  procedure  be  used  to  solve  MGDF  equations. 

The  model  should  be  tested  for  propellants  containing  fuels  other 
than  LP3  and  P13. 

Since  the  MGDF  model  cannot  predict  either  platau  or  mesa 
burning  behaviour  of  propellants,  it  is  suggested  that  further 
modifications  should  be  made  to  enable  the  model  to  predict  these 
experimentally  observed  burning  rate  behaviours. 


■J&sgk] 


113 


REFERENCES 


1.  Summerfield,  M. ,  Sutherland,  G.  S.,  Webb.  M. ,  Taback,  H.J. , 
and  Hall,  K.  P. ,  "Burning  Mechanism  of  Ammonium  Perchlo¬ 
rate  Propellants, "  Solid  Propellant  Rocket  Research,  Progress 
in  Astronautics  and  Rocketry  Series,  Vol.  1,  Academic  Press, 
New  York,  1960,  pp,  141-182. 

2.  Steinz,  J.A, ,  Stang,  P.L.  ,  and  Summerfield.  M. ,  "The 
Burning  Mechanism  of  Ammonium  Perchlorate-Based  Compo¬ 
site  Solid  Propellants,  "  Aerospace  and  Mechanical  Sciences 
Report.  No.  830,  Department  of  Aerospace  and  Mechanical 
Sciences,  Princeton  University,  N.J. ,  1969. 

3.  Williams,  F.  A. ,  Barrere,  M.,  and  Huang,  N.C. ,  "Fundamen¬ 
tal  Aspects  of  Solid  Propellant  Rockets,  "  AGARDOGRAPH  116, 
CIRCA  Publications,  Inc. ,  N.  Y. ,  1969. 

4.  Bakhman,  N.  N. ,  and  Belyaev,  A.  F. ,  "Combustion  of  Hetero¬ 
geneous  Condensed  Systems,  "  Institute  of  Chemical  Physics, 
Moscow  Academy  of  Sciences,  U.  S.S.R.,  Nanka  Publications, 
1967;  also,  RPE  Translation  No.  19  (Editor.  G.  S.  Pearson), 

Nov.  1967. 


114 


5.  Voh  Elbe,  G. ,  King,  M.K.,  McHale,  E.T.,  Macek,  A., 

Friedman,  R. ,  and  Levy,  J.  N. ,  '’Chemical-Kinetics  and 
Physical  Processes  in  Composite  Solid  Propellant  Combustion,” 
Atlantic  Research  Corp. ,  Report  No.  66307,  Jan.  1967.  ) 

i 

6.  Bastress,  E.  K. ,  ’’Modification  of  the  Burning  Rates  of  , 

Ammonium  Perchlorate  Solid  Propellants  by  Partic?e  Size  ; 

Controls,”  Ph.D.  Thesis,  Department  of  Aeronautical 
Engineering,  Princeton  University,  N.J.,  Jan.  1961.  I 

i 

7.  Caveney,  L. ,  Private  Communication,  1972. 

8.  Rumbel,  K.  E. ,  "Poly  (Vinylchloride)  plastisol  Propellants  : 

Pi  opellants  Manufacture,  Hazards,  and  Testing,  "  Advances 
in  Chemistry  Series,  1969,  pp.  36-66. 

9.  Adams,  G.D. ,  Newman,  B.  H. ,  and  Robins,  A.  B. ,  ’’The 
Combustion  of  Propellants  Based  upon  Ammonium  Perchlorate” 
Eighth  Symposium  (International)  on  Combustion,  The  Williams 
and  Wilkins  Co. ,  Baltimore,  1962,  oo.  693-705. 

10.  McAlevy,  R.  F. ,  Lee,  S.  Y. ,  Cole,  R.B. ,  Lastrina,  F.  A. , 
and  Sarnurin,  N.  A. ,  "Investigations  of  the  AP  Composite 
Solid-Propellant  Deflagration  Mechanism  by  Means  of 


Composite  Solid  Propellants, "  Combustion  and  Flame ,  Vol.  12, 
No.  3,  1968,  pp.  201-216. 


|  12.  Culick,  F.  E.C. ,  and  Dehor ity,  G.  L. ,  M  An  Elementary  Calcu- 

i 

i  lation  for  the  Burning  Rate  of  Composite  Solid  Propellants,  " 

j 

j  Combustion  Science  and  Technology,  Vol.  1,  1969. 


|  13.  Wenograd,  J. ,  "Study  of  the  Kinetics  and  Energetics  of 

|  Propellant  Decomposition  Reactions  and  Applications  to  Steady 

i 

< 

State  Combustion  Mechanism,  "  Proceedings  of  the  Third 
ICRPG  Combustion  Conference,  CPIA  Publication  No.  138, 

Vol.  1,  Feb.  1967,  pp.  89-94. 

i 

i 

|  14.  Chaiken,  R.F. ,  and  Anderson,  W.  H. ,  "The Role  of  Binder  in 

|  Composite  Propellant  Combustion,  "  Progress  in  Astronautics 

:  and  Rocketry  Series,  Vol,  1,  Academic  Press,  N.  Y. ,  1960, 

j  pp.  227-249. 


15,  Nachbar,  W, ,  "A  Theoretical  Study  of  the  Burning  of  a  Solid 
Propellant  Sandwich,  "  Solid  Propellant  and  Rocket  Research, 


116 


Progress  in  Aeronautics  and  Rocketry  Series  ,  Vol.  1, 
Academic  Press,  N.  Y. ,  1960,  pp.  207-226. 

16.  Hertzberg,  M. ,  ’’The  Combustion  of  Pure  and  Composite 
Propellants.  The  Expansion  and  Application  of  Laminar  Flame 
Theory  to  Heterogeneous  Solid  Propellants, "  Oxidation  and 
Combustion  Reviews,  Vol.  5,  No.  1,  Cct.  1971. 

17.  Hermance,  C.E.,  " A  Model  of  Composite  Propellant  Combus¬ 
tion  Including  Surface  Heterogeneity  and  Heat  Generation,  " 
AIAA  Journal,  Vol.  4,  No.  9,  1966,  pp.  1620-1637. 

18.  Hermance,  C.  E. ,  "A  Detailed  Model  of  the  Combustion  of 
Composite  Solid  Propellants, "  Proceedings  of  the  2nd  ICRPG/ 
AIAA  Solid  Propulsion  Conference,  CPIA  Publication, 

June  1967,  pp.  89-103. 

19.  Beckstead,  M.W. ,  Derr,  R.L. ,  and  Price,  C.E. ,  "A  Model 
of  Composite  Solid  Propellant  Combustion  Based  on  Multiple 
Flames, "  Lockheed  Propulsion  Co. ,  Report.  No.  TR-69-190, 
1970. 

20.  Beckstead,  M.W. ,  Derr,  R.L. ,  and  Price,  C.E.,  "A  Model 
of  Composite  Solid  Propellant  Combustion  Based  on  Multiple 


117 


Flames, "  AIAA  Journal,  Vol.  8,  No.  12,  1970,  pp.  2200-2207. 

21.  Williams,  F.  A. ,  "Combustion  Theory,  "  Addison  Wesley, 
1965. 

22.  Zelenzik,  F.J. ,  and  Gordon,  S. ,  "A  General  IBM  704  or 
7090  Computer  Program  for  Equilibrium  Compositions,  Rocket 
Performance,  and  Chapman- Jouguet  Detenations,  "  NASA 
TND-1454,  1S66. 

23.  Beckstead,  M.W. ,  Private  Communication,  March  1971. 

24.  Derr,  R.  L. ,  and  Osborn,  J.R, ,  "An  Experimental  Investiga¬ 
tion  of  the  Gaseous  Phase  Reaction  Zone  in  a  Composite  Solid 
Propellant,"  Report,  No.  TM  67-6,  Jet  Propulsion  Center, 
Purdue  University,  Sept.  1967. 


118 


APPENDIX  A  -  FLAME-TEMPERATURE  MEASUREMENTS 

A  modified  Na-D  line  reversal  technique  has  been  used  to  measu¬ 
re  the  flame  temperature  of  the  loose- granule  burner.  The  method 
was  used  previously  by  Sutherland  (1)  and  Derr  (24)  for  similar 
purposes.  The  theory  and  the  derivation  of  equations  have  been 
reported  in  the  cited  references.  Hence,  in  the  following  sections, 
only  the  fabrication  of  loose-granule  burners  suitable  for  temperature 
measurement  is  described  along  with  the  instrumentation  and  error 
evaluation  used. 

This  method  depends  on  measurment  of  the  radiation  intensity 
from  a  flame  in  a  finite  spectral  region,  i.  e. ,  that  of  the  Na-D  lines 
(0.  5890-0.  5896  jui).  Flame  temperatures  were  calculated  from 
variations  produced  by  intermittently  transmitting  light  from  a 
tungsten  lamp  through  the  flame.  Using  Wien's  law  and  Kirchoff  Ts 
law,  tile  flame  temperature  can  be  expressed  as  a  function  of 
1.  brightness  temperature  of  the  tungsten  lamp,  2.  intensity  of  the 
tungsten  lamp,  3.  intensity  of  the  flame  and  4.  intensity  of  the  lamp 
after  being  absorbed  by  the  flame. 

The  following  equation  expresses  this  relationship  : 

J-  =  -L  +  A 

TF  V,u  Ct  Sf 


•  ■*  9 


A-l 


where  T’  is  the  flame  temperature,  T  T  is  the  brightness 
F  Br,  Ij 

temperature  of  the  tujigsten  lamp,  is  the  mean  Na-D  lines  wave 
length,  c2  is  the  physical  constant  in  Wien’s  law  and  D^,  D2,  and  Dg 
are  the  deflections  produced  by  radiation  from  the  lamp,  radiation 
from  the  lamp  after  being  absorbed  by  the  flame,  and  the  ladiation 
from  the  flame,  respectively. 

A. 1  FABRICATION  OF  LOOSE -GRANULE  BURNERS 

Because  the  method  used  is  an  optical  method,  it  requires  that 
the  position  of  the  flame  under  investigation  be  visible  to  a  measuring 
apparatus.  Loose-granule  burners  fabricated  in  lead-foil 
containers  have  been  found  to  be  suitable  for  this  purpose,  since  the 
lead  melted  and  exposed  the  flame  zone  for  temperature  measure¬ 
ments.  The  thickness  of  the  lead  foil  used  was  0. 015". 

A.  2  DESCRIPTION  OF  INSTRUMENTATION 

A  schematic  arrangement  of  the  instrumentation  is  shown  in 
Figure  A-l.  The  apparatus  producing  and  measuring  radiation 
intensities  consisted  of  a  tungsten- filament  lamp,  mechanical 
choppers,  an  interference  filter,  a  photomultiplier,  and  an 
oscillographic  recorder. 
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The  filament  of  the  tungsten  lamp  was  focused  by  the  lens  1  on  a 
# 

hole  in  the  chopper  1.  The  lens  2  focused  the  light  on  the  flame 
and  lens  3‘focused  the  image  of  both  the  flame  and  lamp  filament  on 
chopper  2.  Finally,  the  radiation  was  passed  through  an  interference 
filter  and  focused  on  a  slit  preceding  a  photomultiplier  tube.  The 
resultant  signal  was  amplified  and  recorded  by  an  oscillographic 
recorder. 

When  both  the  choppers  obstructed  the  light  a  zero  reference 
mark  was  obtained;  when  the  holes  in  both  the  choppers  were  open  a 
signal  (D2)  proportional  to  the  intensity  of  the  lamp  radiation  after 
absorption  by  the  flame  was  obtained.  When  f hopper  1  alone 
obstructed  the  light  signal  (Dg)  proportional  to  the  intensity  of  the 
flame  radiation  was  obtained  in  the  absence  of  a  flame. 


A.  3  ERROR  EVALUATION 

Standard  errors  to  be  corrected  for  are  well  documented  by 
Derr  (24) .  All  of  these  have  been  taken  into  account  in  the  present 
investigation.  One  possible  error  which  lias  not  been  considered 
in  the  previous  investigations  (1,  24)  is  that  due  to  signal  noise,  which 
causes  errors  in  the  measurements  of  Dp  Dg,  and  D^.  Using 
Equation  A-l,  the  error  in  T^,,  viz. ,  AT^,  is  found  in  terms  of 
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ADx,  4D2,  and  AD  to  be  : 

4  **  O 


ATf 

X 

r  ad,i 

“  Cl 

i 

i _ 

The  maximum  error  in 

and  AD  are  negative  : 
3 


T  results  when  AD. 
F  1 


is  positive  and  AD 


2 


ATF 


2-^  + 


where  A  D's  represent  the  average  of  AD  ,  AD  ,  andAD  . 

JL  Cl  0 

Typical  magnitudes  of  D  ,  D_,  and  Dc  were  found  to  be  10, 

1  2  3 

8,  and  3  (cm)  respectively,  and  AD  was  observed  to  be  of  the 
order  of  0.  05  (cm).  For  these  values,  the  maximum  error  in  the 
flame  temperature  would  be  of  the  order  of  30°K. 


vmwmmWmm pp 
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APPENDIX  B  -  LISTING  OF  THE  COMPUTER  PROGRAM 
USED  FOR  THE  MGDF-MODEL 
CALCULATIONS 


COMMON  A,  Y/  EF,  RH0/T2/  XLI XR I  /  EPS/  I  END/ Pz  HS/ HI  /  AS>  ESz 

1  FA/FB/D/FC/C/T1/R/T3WI/G/T0/EE/DF/D0X/TAV/ 

2  NJ0B/HAP/CA/CB/CC 
EXTERNAL  FCT 

1  READ<  5/ 1 035NJ0B 
103  FORMAT (I ) 

G0  TO <2/ 3/ 4/5/6/8/9>/NJ0B 

2  READ ( 5/ 105) A/ T2/ TAP/ T0>  P 
READ( 5> 1 06 ) HF/ HAP/ Dr/ D0X 
READ<5/ 107) AS/ ES/ EE/ Y 
READC  5/ 1 08 ) CA/ CB/ CC 
READ( 5/ 109) EPS/ I  END 

GO  TO  1 

3  READ<  5/ 1 05) A/ T2/TAP/T0/ P 
GO  TO  1 

4  READ?  5> 1 06) HF/ HAP/DF  *  D0X 
GO  TO  1 

5  READ<5/ 107)AS/ES/EE/Y 
GO  TO  1 

6  READ(5/ 108)CA/CB/CC 
GO  TO  1 

105  FORMAT ( 5F) 

106  FORMAT (4F) 

107  FORMAT <  4F> 

108  FORMAT <  3F> 

109  FORMAT  <2F) 

8  TAV=(T2+TAP)/2. 

HS=A*HAP+< ! • ~ A) *HF 

RH0= 1 ./<A/1 .95+ C 1 • -A) /DF) 

WRI TEC  6/ 1 10)A/T2/T0/P 

110  FORMAT  <  '  WT/FR= * / F 1 6« 4/  '  T2=’/F10.1/*  T0=*/FlO.l/ 

1  '  P=  ’  /  F  1 0  •  2) 

WRI  TE<  6/  1  1  l>HF/HAP/DF/ASzF.S 

111  FORMAT  C  *  HF=’/F10*2/*  HAP=  •  z  F  1 0 . 2/  ’  DF='/F10.2/ 

1  •  AS= •  /F 1 0»  2/ *  ES= 1  z  F 1  0  »  1  ) 
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WRITE!  6* ! 12)£E> Y/DOX 
112  FORMAT  C  *  EE<~  ' #F20 • 1 0/ '  Y=,/F20.10/  •  D0X= ' > F20. 1 0) 


1 


113 


1  14 


115 


116 


1  17 


12 


«;n 


21 


AB=CA*RH0*TAV**2.75*EXPC22000./C 1.986*TAV))/C 1.95* 
P**2.*D0X**2. ) 

AF=A8/< l.+AB) 

AC=CA 

AD=CB 

WRITEC6/  1  13)AC/AD>AB 
FORMAT ( '  CA=  *  * E20 .10/' 

WRITEC6#  1  1/OAF 
FORMAT C  '  WT  FR  FUEL**, 


CB= ’ *  E20. 1  Of '  Ms  S  E20. 10) 


20 


F40. 10) 

T0NE= 1384.687+32252. 1 2* AF-23 5047. 4*AF**2. + 7021 79* 
AF**3. 

WRITE  C  6/ 115) TONE 
FORMAT C  *  TOt\'E=  * /F20 » 5) 

HI  1  =  0*  3* ( TONE -TO) +HAP 
HI =-A*Ml 1 
WRITE <6# 1 1 6)HI 
FORMAT ! ’  His ' >F20. 5) 

TOW I =CB*TAV*EXPC 22000*/ ( 1 *986*TAV) >/P 
WRITE C6/ 1 17) TOW I 
FORMAT <  '  T0WI  =  SE3O»  10) 

E=EE 

EF=E*T2*EXP!20000«./C2»*1  *986*T2))/P 
GO  TO  1 

READ <5# 12>XLI/XRI 
FORMAT (2F) 

CALL  RTMI  CX/F/FCT/XLI/XRI/EPS/IEND/IER) 

IF <IER.NE*O)G0  TO  101 

WRITE!6/50)R 

FORMAT ! /////£ 1 6 18  ? 

T0W2sP*3O**!EF+G>**2./C82.O5*T2> 

EL2s82.O5*RH0*R*T2*TOW2/C3O.*P) 

EL1=82.O5*RH0*R*T1*T0WI/C3O.*P) 

OT=C 

QT1=R*RH0*OT 

QII=.0002*!T2-T1 )/<EL2+RH0*R> 

QI=QT-OII 
01 1=R*RH0*OI 
01  IF=0I I*R*RH0 
TOW  1 sTOW I 

W  R I T  E ( 6  >  2 1 ) TO W 1 /  T  0  W  2 / EL \» ELS 
FORMAT! ////E:i  6.8^  El  6. 8^  El  6.8/ El  6*  8) 

WRI  TEC  6/  20)  QT/QT 1 

FORMAT! *  0T= ' >  E20 .6/ ’  QTF=‘/E20.6) 


# 


see  note  on  page  125 
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WRITEC6»22)QI?QI  ! 

22  F0RMATC  *  QI-,?E20.6?,  GIF=*?E20.6) 

WRITE  C6?23)GII?QIIF 

23  F0RMATC*  QI 1  =  ' ? £  1 0? 4?  ’  01 1  F=  *  »  E 1 0.  4> 

ET=P/CC 

EV=P***6666666 
WRI TEC  6?  300) ET?  £V 

300  FORMAT C  *  P/R=  ' ? F 1 0*  5?  *  P2/3s= v ?F 1 0»  5> 

GO  T0  1 

101  WRITE (6* 1 0  ) I ER 

100  F0RMAT C  *  ERROR  C0DE=SI3> 

G0  T3  9 
END 

FUNCTION  FCTCX) 

COMMON  A?  Y>  EF?  RM0?  T2>  XLI  >  XRI  *  EPS>  I  END?  P?  HS?  HI  ?  AS>  F.S? 

1  FA?  FD?D?FC?C?T1?R?T0WI?G?T0/>EE?DF?D0X?  TAV? 

2  NJ0B?  HAP?  CA?  CB?  CC 
R-AS*EXPC-ES/C 1 .986+X) > 

C2=0 

C«.3*CX-T0)+HS 

D=<0002*CT2-X) 

G= C Y*T2**. 83333) /CX**. 87 S*P**. 3333333) 

WRITEC  6?2>X?  R 

2  FORMAT <// *  T5= ' ?  F 1 2 • 4?  *  R=’?E16*6> 

FA=R**2. 

FB=1  .-CFA*T0WI  *82. O5*RH0**2*C? / C P*0 ?  0002*30*) 
FC=-.3*T0+HS+HI*C2 

FCT=F A-  .  0002*  CT  P-X/FB )  /  C  C  EF  +  G  )  **2 .  *RH0**2  . 

1  *C .3*X/fB+FC) > 

RR=SQRTCFA-FCT) 

T 1 =X/FB 

c:  \  *r  *  no 
tt  m  H-  v  u /  J  /  i  i  ;  no 

5  FORMAT C  *  T 1= ’ ?F 12 .4? '  RR=’?E16«8) 

CC=R 

RETURN 

END 


Note  :  * 

T.  (TONE)  given  in  this  program  is  for  AP/P13 
propellant  (valid  in  0.  9<  w  <  1. 0  ) 

* 

T  (TONE)  for  AP/LP3  propellant  is  given  by 
TONE  =  1398.  21  +  13891.  07  *  AF  -  33482. 12*AF**2.  0 


